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^ \ Abstract 
> ■ 

\^ • The fundamental theory of nature may allow a large landscape of vacua. Even if the 

'nJ" ■ theory contains a unified gauge symmetry, the 22 flavor parameters of the Standard Model, 

including neutrino masses, may be largely determined by the statistics of this landscape, 
and not by any symmetry. Then the measured values of the flavor parameters do not 
^ ■ lead to any fundamental symmetries, but are statistical accidents; their precise values do 

! not provide any insights into the fundamental theory, rather the overall pattern of fla- 

vor reflects the underlying landscape. We investigate whether random selection from the 
statistics of a simple landscape can explain the broad patterns of quark, charged lepton, 
^ , and neutrino masses and mixings. We propose Gaussian landscapes as simplified models 

\ of landscapes where Yukawa couplings result from overlap integrals of zero-mode wave- 

functions in higher-dimensional supersymmetric gauge theories. In terms of just five free 
parameters, such landscapes can account for all gross features of flavor, including: the 
hierarchy of quark and charged lepton masses; small quark mixing angles in the basis with 
quarks arranged according to mass, with 13 mixing less than 12 and 23 mixing; very light 
Majorana neutrino masses, with the solar to atmospheric neutrino mass ratio consistent 
with data; distributions for leptonic mixings sin 2^12 and sin 2^23 that are peaked at large 
values, while the distribution for sin 2^13 is peaked at low values; and order unity CP vio- 
lating phases in both the quark and lepton sectors. While the statistical distributions for 
flavor parameters are broad, the distributions are robust to changes in the geometry of the 
extra dimensions. Constraining the distributions by loose cuts about observed values leads 
to narrower distributions for neutrino measurements of ^13, CP violation, and neutrinoless 
double beta decay. 
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Table 1: The 17 measured flavor parameters and a limit on a mixing angle in the lepton sector. 
All data comes from the pdgLive feature from the Particle Data Group [1]. For comparison to 
distributions provided throughout this paper, the quark and charged lepton Yukawa eigenvalues 
and the CKM matrix elements have been run up to the (reduced) Planck scale Mp = 2.4 x 10^^ 
GeV, assuming no physics beyond the Standard Model enters up to this scale. On the other 
hand, the RG scaling effects are not taken into account in the neutrino sector. 
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1 Introduction 

The Standard Models] describes all laboratory data with 28 free parameters. Of these, 22 arise 
from the Yukawa matrices A"''''*^ and the coupling matrix C that appear in the flavor interactions 

C 

l^fiav = KjUiqj h + XtjdkQj h* + XliCji h* + -j^kljhh , (1) 

where q, I {u, d, e) are the left (right) handed quark and lepton fields. Of these 22 flavor param- 
eters, 17 are measured, with varying levels of accuracy. The remaining five parameters are all 
in the neutrino sector and two of these have upper limits. Given several decades of continued 
progress on improving the accuracy of the experimental measurements, the most striking fact is 
that there is nothing approaching a standard theory of the origin of these parameters. Not only 
are we ignorant about the overall picture of flavor, we do not have a convincing explanation for 
the value of any of the 17 measured parameters. 

A complete theory of flavor would provide answers to three very different questions: 

• What is the origin of the fermion quantum numbers and why are there three generations? 

• What determines the qualitative pattern of the quark and lepton mass matrices? For 
example, why do the charged fermion masses and mixings have a hierarchical pattern, 
while in the neutrino sector there are large mixing angles? 

• What determines the precise values of the 22 flavor parameters? 

The couplings A"'*^'*^ and C are symmetry breaking parameters of the flavor symmetry group 
U (3)^, where one U (3) factor acts on each of q, I, u, d, e. The dominant approach to constructing 
theories of flavor is to use symmetries to reduce the number of free parameters, n. An underlying 
flavor group Gf G U{3)^ and a more unified gauge symmetry both limit n, leading to precise 
predictions if n < Uobs, the number of observables. A hierarchy of symmetry breaking scales can 
lead to small dimensionless parameters that explain qualitative features of the mass matrices 
[2], and there are many realizations with Gf both Abelian and non-Abelian. Still, it is striking 
that the progress along these lines is limited, even as more precise data have become available. 
Perhaps this is a sign that a completely new approach is needed. 

The cosmo logical dark energy [3] apparently has little to do with flavor. However, it may 
be the first evidence for a huge landscape of vacua, with the observed value of the cosmological 
constant resulting from environmental selection for large scale structure [1]. If we take string 
theory to be a theory of a landscape rather than a theory of a single vacuum, then what are 
the implications for flavor? The enormous number of vacua, required for a sufficiently fine scan 
of the cosmological constant, results from the large number of ways background fluxes can be 
arranged on non-trivial compact manifolds of extra spatial dimensions [S]. The cosmological 
constant, however, is not the only parameter that scans over such a landscape. The Yukawa 
couplings may also vary from one universe to another. 

^We everywhere refer to the "Standard Model" as the theory including the dimension-five operators of ([T]). 
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In any landscape of vacua, the Yukawa couplings have a relevant dependence on n = ns + nF 
parameters, where of the parameters scan and np oi them are fixed. If < Uobs then data can 
determine the subset of vacua in which we happen to live and nobs~n predictions can be made. A 
certain simple model with a single extra dimension has Uobs — n = 1, giving a single prediction [B]. 
If n > Hobs, however, then no such precise predictions are possible. Some flavor parameters — such 
as the electron, up and down masses — may be determined by environmental effects [7j. The top 
mass may also be determined along this line [8], while selecting for leptogenesis may be a key 
to understanding some parameters in the lepton sector. However, most of the flavor parameters 
do not seem to be strongly selected. Therefore it could be that most flavor parameters have 
values that simply reflect some underlying probability distribution over an enormous number of 
solutions to the fundamental theory. The precise value of any such observable is accidental and 
not fundamental, since any other nearby value would be just as probable. Although it may be 
unappealing that there is nothing fundamental or beautiful relating the flavor parameters in our 
universe, so far this possibility cannot be dismissed, especially given that no such relations have 
been found. 

The large mixing angles observed in atmospheric and solar neutrino oscillations inspired the 
idea that the relevant Yukawa couplings of the neutrino sector are governed by randomness rather 
than by flavor symmetries [H] . By introducing a simple probability distribution for the elements 
in the Dirac and Majorana neutrino mass matrices, probability distributions for the neutrino 
observables were generated and found to agree well with data. However, a complete landscape 
approach to the theory of flavor must explain many other things; there are rich structures among 
the observed flavor parameters. Reference [10] took a bold step in this direction, studying 
masses and mixings in both the quark and lepton sectors by introducing a simple probability 
distribution for the Yukawa couplings. We begin in section [2] by analyzing the toy landscape 
introduced in [T^. While the main purpose of this section is to clarify how hierarchy among 
quark masses is generated in this model and to point out its limitations, we also introduce an 
approximate analytical study of the probability distribution of flavor observables. This turns 
out to be a useful warm-up exercise for the subsequent sections. 

Ideally, a statistical approach to the theory of flavor would be both phenomenologically suc- 
cessful and theoretically well-motivated. However to determine the probability distributions for 
flavor parameters from a purely top-down calculation is not an easy task; meanwhile the number 
of observables is limited, and such low-energy information does not constrain the underlying 
probability distributions in a purely bottom-up approach. We consider as a practical strategy 
to use string theory as a guide in deducing the form of the probability distributions, while using 
known experimental constraints to try to find a phenomenologically successful model within 
this restricted set of possibilities. In a given compactification of string theory, the Yukawa cou- 
plings are determined by an overlap integration over the extra spatial dimensions. Scanning 
moduli parameters of the gauge field configuration on the extra dimensions corresponds to scan- 
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ning zero-mode wavefunctions, and hence this corresponds to scanning Yukawa couphngs. If 
the wavefunctions of the quarks, leptons, and the Higgs are peaked at different locations, then 
the overlap integral can lead to small Yukawa matrix elements. Hence localization in extra 
dimensions is an alternative to symmetries for generating fermion mass hierarchies [TT]. 

In sections [3] through [6] we introduce simple "Gaussian landscapes" and study the resulting 
distributions for quark and lepton masses and mixingsjfl These landscapes have features expected 
from certain string landscapes, for which they can be viewed as simplified or toy models. The 
key feature of a Gaussian landscape is that all quark, lepton and Higgs fields have zero-mode 
wavefunctions with Gaussian profiles in the extra dimensions, and that the centers of these 
profiles all scan independently with fiat probability distributions over the volume of the extra 
dimensions. For simplicity, the number of free parameters used to describe the geometry of the 
extra dimensions and the widths of the Gaussian profiles is kept to a minimum. We find the 
observed quark and lepton masses and mixings can be typical in such simple landscapes. 

In section [3] we introduce a Gaussian landscape for quarks on a circle S^, with all Gaussian 
profiles having the same width. Numerical probability distributions for the nine CP-conserving 
flavor observables are provided, and a qualitative semi-analytic description of these distributions 
is derived. The results are compared with those that result from introducing approximate flavor 
symmetries, with some similarities and some differences emerging. Finally, the effects of possible 
environmental selection on the top quark mass is studied. The large number of flavor parameters 
in the Standard Model allows for a reasonably significant evaluation of goodness-of-fit between 
a Gaussian landscape and the observed flavor structure. This is described in section HI with the 

Gaussian landscape used for illustration. The effects on the quark sector from adding more 
dimensions to the Gaussian landscape are examined in section [5l together with a preliminary 
study of the effects of geometry. 

In section [B] the Gaussian landscape is extended to include the lepton sector. Motivated 
by the expectation of supersymmetry in the higher dimensional theory, Yukawa couplings for 
both Dirac and Majorana neutrino masses are introduced in the Gaussian landscape, with light 
neutrinos resulting from the usual seesaw mechanism in four dimensions. Larger mixing angles 
in the lepton sector might arise from assigning the appropriate fermions larger Gaussian widths, 
but these are found to largely cancel between charged and neutral sectors. However, this cancel- 
lation is prevented by introducing CP-violating phases in the Gaussian profiles. This suggests a 
connection in Gaussian landscapes between large CP violation and large leptonic mixing angles. 
In addition to providing numerical distributions for the observed lepton flavor parameters, prob- 
ability distributions are predicted for the leptonic mixing angle ^13, the CP phase in neutrino 
oscillations, and the Majorana mass relevant for neutrinoless double beta decay, illustrating how 
these toy models for the string landscape can connect to future experimental data. Finally, the 
extent to which these Gaussian landscapes arise from unified super symmetric field theories in 

brief summary of some of the major results of this work can be found in [12j . 
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higher dimensions is examined in section [3 Such theories have huge symmetries — ultimately all 
entries of the Yukawa matrices arise from a unified supersymmetric gauge coupling — so that the 
complicated pattern of observed masses and mixings can arise from a very simple mechanism: 
the scanning of the centers of the Gaussian profiles. 



2 Prelude: Hierarchy without Flavor Symmetry 

In previous work it has been suggested that the components of the Yukawa matrices A"'*^'^ and the 
coupling matrix C are selected randomly and independently of each other [HI HO]- For example, 
in Neutrino Anarchy [9] one finds that the large mixing angles underlying neutrino oscillation 
are typical of the lepton interactions that arise when each element of the matrices A*^ and C is 
independently selected from the simple probability distributions 

dPiY) dPiC) 

= const., = const., (2) 

or from distributions such as 

'^^}\] = const, for A^j^ < A^ < A;;,^^ , = const, for Cmin <C < C^ax , (3) 

a m A"^ amC 

where X'^ij^^^.x C*min,max are of order unity (for more details see P [13]). Such an absolute 
anarchy of lepton couplings tends to result in comparable mass eigenvalues. On the other 
hand, [TU] introduced a power-law probability distribution for the Yukawa matrix elements, 
dP{X)/d\ oc A""^ for Amin < A < Amax and dP{X)/dX = otherwise. By assuming Amin <^ Xma.x 
and choosing 6 appropriately, quark Yukawa matrices with each matrix element following such 
a distribution can roughly accommodate the hierarchical pattern of quark mass eigenvalues |10j . 
According to [1^, 6 = 1.16 provides the best fit to the quark sector. 

In this paper we propose a significant modification to these ideas. However as a prelude to 
discussing our proposal it is worthwhile to first study the model of [10] in greater detail. To 
simplify the analysis we specialize to the particular case 6=1, such that 



dP{X) 
d logio 




niiax 



for A > A„ 

for Amin < A < Amax • (4) 
for A < Amin 



Henceforth we refer to this distribution as a scale-invariant distribution. In this section we 
also restrict attention to the quark sector. Note that CP-violating phases are not introduced 
in this landscape because all of the matrix elements are real-valued. A pair of 3 x 3 Yukawa 
matrices A" and A'^ is generated by choosing each of the 18 matrix elements randomly according 
to the distribution (jl]). The Yukawa matrices are then diagonalized using separate left- and 
right-handed unitary transformations, and the quark masses and mixings are calculated from 
the eigenvalues of the diagonalized Yukawa matrices and the resulting CKM matrix. The above 
process is then repeated to generate an ensemble of these nine observables. 
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Figure 1: Distributions of the three eigenvalues of Yukawa matrices whose each element follows 
the distribution (jlj). From left to right the three panels correspond to the smallest, middle and 
largest eigenvalues. The sum of all three distributions reproduces Figure 9a of |TU]. We used 
log^Q Ajnin = —9 and log^Q Amax = 0.1 for this simulation. 



2.1 The Distribution of Mass Eigenvalues 

We first study the distributions of the quark masses. Since the up- and down-type Yukawa 
matrices are generated independently of each other and in exactly the same way, the distributions 
of masses for these sectors are exactly the same. Thus we only need to study one of the two 
sectors. Results of a numerical study are shown in Figure [H where we have chosen \min/ Amax = 
10~^-^. With some approximations we can understand the shapes of these distributions. Let A3 
denote the largest element of the 3x3 matrix A. Meanwhile, the largest element of the 2x2 
sub-matrix of A that excludes A3 is denoted Ag. For example, in the matrix 

(5) 





( A23 


A2 A21 




Al3 


A 12 A'l 




V ^3 


A32 A31 



A3 is the largest among the nine entries and Ag the largest in the upper right 2x2 sub-matrix. 
Given this characterization, the probability distribution for the variables A'^^ 2 3 ^^'^ ^ij is 

d-P{x[,2,3^ Xij) = 36 9(x3 - X2)Q{x'^ - a;32)e(a;3 - X23)&{x'3 - X3i)Q{x'^ - Xi3)6(x2 - x[) 

X0(X2 — X2l)0(x2 — 0:12) dx'idx'2dx'3dXi2dX2ldXisdXsidx 23(1X32 , (6) 

where is the step function. The factor of 36 comes from the nine possible locations for A3 
times the four possible locations for Ag. In addition we have introduced the notation 



^/ _ lll(A'j/Amin) X- — (J'j 

ln(Amax/ Amin) ln(Aniax/ An ^ 



The largest mass eigenvalue of is approximately A3; this approximation is poor if one of 
A325 A31, A23 and Ai3 is almost as large A3, but this is unlikely if ln(Amax/Amin) is large. We call 
this largest eigenvalue A3 and define X3 analogously to Xg. The probability distribution of X3 is 
therefore approximated by integrating out from all the variables except x'^: 

rfP(x3) ~ dX3 . (8) 
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Figure 2: The approximate distributions of the three eigenvalues given in flH HTOjl . and those of 
mixing angles given in (fT5| [T6l) . We have used the same Amax and Amin as are used in Figured! 



Meanwhile, we approximate the middle eigenvalue A2 by Ag and the smallest eigenvalue Ai by 
X[. This approximation is poor when the seesaw contributions (Aj3A3j)/A3 and (Ai2A2i)/A2 are 
larger than Aj and A'^, respectively. Thus we do not expect this approximation to be reliable for 
small values of X2 and Xi. Nevertheless, integrating out all of other variables we find 

36 

dP{x2) = —xl{l-xl)dx2, (9) 
5 

dPix,) = f (10) 

The probability distributions (P 4T0I) are shown in Figure O Remarkably, they capture the gross 
features of the numerical results in Figure [H Therefore we use these distributions to examine 
the qualitative aspects of the mass distributions that follow from the landscape (jl]). 

The average (xj) and the standard deviation of the three Yukawa eigenvalues (both on a 



logarithmic scale) can be calculated from the distributions fl8]- [T0|l : 

(X3) = 0.90, da = 0.09, X3 ~ [0.81 - 0.99] , (11) 

(X2) = 0.72, ^2 = 0.16, X2 ~ [0.56 - 0.88] , (12) 

(xi) = 0.36, ai = 0.22, xi ~ [0.14 - 0.58] . (13) 



The three eigenvalues are on average well-separated and they overlap with neighboring eigenval- 
ues only slightly at one standard deviation. Even this slight overlap between the distributions 
is misleading. Recall that by definition for any particular set of Yukawa matrices we have 
X3 > X2 > x'l- Thus the combined distribution for the eigenvalues is not a naive product of 
(IH HTOl) but is given by integrating the other six variables out of ([6]). This gives 

dP{xi, X2, X3) ~ 36 x'^xl 9(x3 — X2)6(x2 — xi) dxidx2dx3 . (14) 

Thus it happens that in the subset of cases where x^ is small, the distribution of X2 is pushed 
to even smaller values. Hence the three eigenvalues tend to be well-separated even in the log- 
arithmic scale, and only rarely are adjacent eigenvalues comparable. Note that none of this 



8 



. 1 0.2 



logio(2ei2/vr) 



. 1 0.2 0.3 0.4 



logio(2^23/7r) 



JJ 



25000 
20000 



30000 
25000 
20000 
15000 



0.4 0.6 0.8 



sin 9 



13 




logio sin 6*13 



Figure 3: Distributions of the three mixing angles of the CKM matrix that results from the 
distribution (jlj). The upper left panel is similar to Figure 11 of [TU]; however in [TU] 5 = 1.16 is 
used and the distribution is displayed only within the range < 6ij < 7r/4. 



depends on the choice of Amin and Amax- In short, hierarchical structure (Yukawa eigenvalues 
well-separated in logarithmic scale) is generated statistically in a landscape where each matrix 
element independently follows the scale-invariant distribution (j4]). Whether the hierarchy is 
large or small is determined by whether logiQ(Amax/Amm) is large or small. 

2.2 Pairing Structure in Electroweak Interactions 

Let us now study the mixing angles. Figure [3] shows the distributions of mixing angles in the 
quark sector that result from a numerical simulation where each element of both the up-type 
and down-type Yukawa matrices is assumed to follow the distribution (j4]) independently^ The 
probabihty distribution functions of the mixing angles are shown against the axes of d9i2, (1023 
and d{sm 613) because the invariant measur^ of S0(3) mixing matrices is d9i2 A d923 A d{sm 613). 
Since in this model the Yukawa couplings are all real- and positive-valued, the CKM quark 
mixing matrices are SO (3) matrices. 

The prominent feature of these distributions is the twin peaks at 6ij = and 6ij = 7r/2 for all 

^ Using ([5 HTU)) we find these distribution functions to be given by the approximate analytic form: 

dP{t) ~ — {10- im + 6t^ -t"^) {5- 8t^ + 31"^) dt, where t= , , , (15) 

50 ln(Arain/Amax) 

dP{t) Ii8-9t + f){l-f)dt^ whereto Z.^^^.^^-^^ . (16) 

"^The importance of the invariant measure is emphasized in [13] , 
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three mixing angles. This feature is straightforward to understand. Suppose that the randomly 
generated Yukawa matrices are of the form 

/ * A- . \ / . * A'^\ 

~ * • A'H , Xtj-\ * * ], (17) 

V A- * * ; V A'? * . / 

where the *'s are assumed to be less than A'g or A'3, and the w's less than A'2 or X"^- Ignoring 
the seesaw contributions to eigenvalues, i.e. when (* * /A3) <^ Ag and (• • /Aj) <^ A'^, we find 
that the three left-handed quark doublets qj (j = 1,2,3) are approximately 

qi = ith, di), q2 = (cl, &l), qs = (ul, Sl), (18) 

where t^, cl and ul stand for left-handed components of the heaviest, middle and lightest 
mass eigenstates of up-type quarks. The down-type mass eigenstates 6^, sl and are defined 
similarly. The CKM matrix for the up-type and down-type Yukawa matrices (fTTll is roughly 



VcKM ~ 1^ 1 j , (19) 

corresponding to 612 ~ ^23 ~ 7r/2 and ^13 ~ 0. This explains one of the 2^ combinations of 
peaks in Figure [31 When the \%^2i found in different entries of the 3x3 Yukawa matrices, 
other peak combinations are obtained. Continuous distributions connecting 9ij ~ to 6'^^ ~ 7r/2 
originate from the seesaw contributions that we have ignored. Thus the distributions become 
more and more localized around the peaks as Amax/Amin is increased and seesaw contributions 
become less important. Therefore the small mixing angles of the observed CKM matrix are not 
atypical of distributions with large values of Amax/Amin- 

How then do we interpret the peaks at 7r/2? The flavor structure of the quark sector of 
the Standard Model is characterized by three general features. On the one hand, the quark 
masses have a hierarchical structure, and this is successfully reproduced. In addition, the W- 
boson current approximately connects three distinct pairs of quarks — we refer to this as "pairing 
structure," and it is also found in this landscape. Finally, the VT-boson current connects pairs 
such that the lightest up-type quark is approximately paired with the lightest down-type quark, 
the middle up-type quark is approximately paired with the middle down-type, etc. We refer to 
this as the "generation structure" of the Standard Model. Mixing angles near 7r/2 maximally 
violate this generation structure; for example the set of angles 612 ~ ^23 ~ 7r/2 and 6*13 ~ of 
fll9p corresponds to the quark pairings in fll8p . There are 3! combinations in forming three pairs, 
only one of which has what we call the generation structurejfl 



^ When sin 013 ~ 1, 612 ~ ^23 ^ and 612 ^ 623 ^ tt/2 result in the same quark combinations. Likewise, 
(012,6*23) ~ (0i'''/2) and (0i2,023) ~ ("""/SiO) are the same. Therefore there are only 3! physically different 
combinations, even though there are 2^ different ways to pick three peaks from Figured 
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2.3 Problems 



Is the existence of mixing angle peaks about 9ij = 7r/2 really a problem? The landscape that we 
have discussed so far may reproduce the generation structure of mass eigenstates in the H^-boson 
current (when 612 ~ 6*23 ~ 6*13 ~ 0), but more often it does not. Although it might be argued 
that this is just a 1/3! coincidence problem, it is still difficult to accept that the generation 
structure of flavor is not revealing something important about the underlying theory. It is also 
tempting to try read something deeper from the observed hierarchy between Vub ~ ^13 ~ 4 x 10~^ 
and the Cabibbo angle 9i2 ~ 0.2. In this landscape this hierarchy is just a random statistical 
fluctuation (c.f. Figures [2] and [3]) , and it does not appear this shortcoming can be overcome by 
a more ideal choice of Amax/Amin (c.f. Figure [1]). 

An even bigger problem is to understand how the probability distribution (jl]) arises or what 
is the correct distribution to replace it. That is, although the phenomenology of the landscape 
that we have considered may be deemed acceptable, we do not have a solid theoretical ground 
upon which to base it|§ It is also obscure how the phenomenology of a lepton sector with large 
mixing angles and that of quark sector with small mixing angles can be accommodated within a 
single theoretical framework. In the remaining sections of this paper we analyze some landscape 
models that successfully reproduce the phenomenology of hierarchy, pairing, and generation 
structure, while making progress on each of the four problems described above. 

3 A Gaussian Landscape for Quarks in One Extra Di- 
mension 

The landscape discussed in section [2] assumes that all 18 elements of A" and A*^ are scanned 
independently. Yet without any correlation between these two Yukawa matrices, generation 
structure will never be obtained. For example, in order to ensure that the heaviest up-type 
quark ti is contained in the same SU(2)l doublet as the heaviest down-type quark 6^, we 
require the following. When an {i, j) element of the up- type Yukawa coupling X^^ Ui qj h is large, 
at least one of the three down-type Yukawa couplings Xfj dk qj h* {k = 1, 2, 3) involving the same 
quark doublet qj should be large. A landscape of vacua must realize such a correlation between 
the up-type and down-type Yukawa matrices in order to explain the generation structure. 

Perhaps one of the simplest ideas to introduce such a correlation is to introduce an extra 
dimension. If a large Yukawa coupling of A^- Ui qj h is due to a substantial overlap of the wave- 

^Such an attempt is made in |14j . However, we consider that correlation among Yukawa couplings and the 
number of extra dimensions are crucial ingredients in understanding flavor physics, and these are missing in [14j . 
The intersecting D6-D6 system mentioned in [T3] is dual to a T^-fibered compactification of Heterotic string 
theory and is simulated by the D = 3 Gaussian landscape models of this article. Although [Mj guesses that 
the scale invariant (or nearly scale invariant power-law) distribution might arise from the intersecting D-brane 
systems, we conclude otherwise. Specifically, we find that the scale-invariant distribution ([4]) is derived from the 
D = 1 Gaussian landscape model, but not a D = 3 model. We also find that correlation among various elements 
of the Yukawa matrices is crucial to understanding flavor structure and the origin of generations. 
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functions of qj and h, then the down-type Yukawa couphngs involving the same qj tend to be 
larger because of the overlap of qj and h. At the same time, for localized wavefunctions the 
overlap of some triplets of u, q and h can be very small, and so there is hope to explain the 
hierarchically small Yukawa couplings necessary to account for light quarks. 

In this section we present a simple toy landscape based on Gaussian wavefunctions spanning 
a circular extra dimension. Although a single extra dimension is introduced for simplicity, this 
model captures the essence of what one expects more generally from such "Gaussian landscapes" 
based on multiple extra-dimensional field theories. Through numerical simulation and an ap- 
proximate analytical analysis, we find that the hierarchy, pairing, and generation structures of 
quarks is obtained statistically in this landscape. No flavor symmetry is needed. 

3.1 Emergence of Scale- Invariant Distributions 

We introduce a single extra dimension with the simplest geometry: S^. The wavefunctions for 
all of the quarks and the Higgs boson are assumed to be Gaussian with a common width d, and 
centered at arbitrary points on S^: 

1 (y-yp)^ 

Here y is the coordinate of and M5 is the cut-off scale of the effective field theory in 4 -|- 1 
dimensions. This wavefunction is normalized so that 

M5 / dyip\y) = l, (21) 







where L is the circumference of S^. The wavefunction ( l20i) should be made periodic on S^, 
while maintaining the normalization in (^T^. Yet as long as the width of the Gaussian profile d 
is parametrically smaller than the circumference L, the wavefunction is almost Gaussian. One 
should examine whether Gaussian wavefunctions arise as solutions to equations of motion of field 
theories in extra dimensions, but we defer this theoretical study to section [3, and first study 
whether the assumption of Gaussian wavefunctions on extra dimensions leads to a successful 
explanation of the physics of quark and lepton masses and mixing angles. 

We calculate the up-type and down-type Yukawa matrices with the overlap integrals 



A" = gM, dy^f{y;a,)^]{y;bj)^\y;y'^), 

Af,. = gM, [ dy4{y-Ck)^]{y,h,)^\y-y'^), (22) 

where g is an overall constant Jl| ^'j{y), ^iiy), ^tiv) V^iv) ^-^^ wavefunctions of left-handed 
quark doublets q^ (j = 1, 2, 3), right-handed up-quarks Ui {i = 1, 2, 3), right-handed down-quarks 

^We will see in section [7] that these interactions may originate from super Yang-Mills interactions on a higher 
dimensional spacetime, and then g derives from the Yang-Mills coupling constant. Despite this origin, g can be 
different in different sectors in effective descriptions such as those using just one extra dimension. For simplicity 
wc set the constant g to be the same for both the up and down sectors. 
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dk [k = 1, 2, 3) and of the Higgs boson, respectively, all of the Gaussian form fl20|) FI The center 
coordinates of these wavefunctions are 6j, Oj, Ck and y^, respectively. The matrices A"''^ are 
real, so that CP is conserved in this toy landscape. Complex Gaussian wavefunctions will be 
studied in section Ej where the effects of phases on the distributions of quark masses and mixing 
angles is found to be small. We assume that the center coordinates 6j, a^, and are scanned 
freely and independently from one another on . Because of the translational symmetry of 
S"^, only the relative difference between these center coordinates affects observables. Thus the 
effective number of scanning parameters is ns = 9. On the other hand, there are nine observables 
determined from the Yukawa matrices in the quark sector: six mass eigenvalues and three mixing 
angles. Thus the scanning parameters cover the space of observables and no precise prediction 
among the observables is available. However, since this Gaussian landscape covers the space of 
observables, our vacuum is unlikely to be missed in this ensemble. 

The other model parameters, namely the width d, circumference L, cut-off scale M5, and 
coupling g, are treated as fixed. This treatment is quite arbitrary; among the myriad of other 
possibilities are to scan some or all of these parameters, to allow the up-sector and the down- 
sector to have different values of g, to choose different widths d for different wavefunctions, 
etc. An extreme version of the landscape would allow everything to scan, leaving no fixed 
parameters to be input by hand. Our treatment — namely four fixed parameters — is equivalent 
to shcing a specific subset out of a possibly much larger landscape. This allows us to identify a 
phenomenologically successful subset of the landscape and at the same time more easily explore 
its properties. It is a separate question whether this subset is typical of the full landscape, 
or if it is highly weighted because of cosmological evolution and/or environmental selection. 
Indeed, it is not implausible that the distributions of these four fixed parameters could be 
sharply peaked in the full landscape, since some toy landscapes predict Gaussian distributions 
for some parameters and cosmological evolution can provide exponentially steep weight factors. 
We consider that practical progress can be made by dividing the full problem into simpler parts. 

Note that out of these four fixed parameters there are only two independent combinations 
that affect the Yukawa matrices. First, only the three dimensionless combinations, g, M^d, and 
M5L, are physical; the value of M5 simply sets the scale for dimensionful parameters. Second, 
even as we change the "volume" of the extra dimension M5L, the Yukawa couplings remain the 
same if the width parameter M^d and the coupling g are scaled so that the ratio M^d/M^L = d/L 
and the effective coupling, 

9cs = , ^ , (23) 



*In sectioii[7]we study higher dimensional supersymmetric field theories, where at short distances the up-type 
and down-type Yukawa couplings involve difi^erent Higgs doublets, /ii,2- Our results are independent of whether 
or not supersymmetry survives to the weak scale. If a supersymmetry survives compactification, we assume that 
the hi^2 zero modes have the same localization. If no supersymmetry survives then we assume a single Higgs 
zero mode with a Gaussian profile. There is a "tan/3" factor between the up and down sectors that we have set 
to unity. With weak scale supersymmetry this factor arises from the ratio of Higgs vacuum expectation values, 
while with high scale supersymmetry breaking it arises from the composition of the light Higgs boson. 
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remain the same. Therefore the Yukawa couphngs are effectively controlled by only these two 
parameters, d/L and Qcq. (Note however that the volume M5L does affect the low-energy value 
of Newton's constant.) 

The Yukawa couplings, given by the overlap integrals (1221) . can be expressed more explicitly 
in terms of the underlying parameters in a restricted region of the parameter space. Suppose 
that d/L ^ 1. Then the compactness of is not important in the calculation of the Yukawa 
couplings, as long as the center coordinates of quarks, (or Ck) and 6j, are close to that of the 
Higgs boson (which, by translational invariance, we set as the origin of the coordinate y). For 
such a vacuum, the Yukawa couplings are given by 



Let us compare this result to the form for the Yukawa couplings that results from approximate 
Abelian Flavor Symmetries (AFS). In the most general AFS scheme there is a small symmetry 
breaking factor associated with each quark field, e'^'"''^, which leads to Yukawa matrix elements 

Al^ = 4^ , = g'uj 4^ , (25) 

where the gij and g'^.^ are all of order unity. A mass hierarchy among the generations is realized 
by imposing ei in the left, right, or both sectors. Models with fewer parameters can 

be constructed and then the symmetry breaking parameters are not all independent; consider 
for example a single Abelian symmetry with a symmetry breaking parameter e that appears in 
different entries with different powers due to generation dependent charges. Generation charges 
(0,2,3) then give €3 : 62 : ei = 1 : : e^. No matter how the model is arranged, the mass 
hierarchy arises because the first generation feels much less flavor symmetry breaking than the 
third. Note that AFS theories are very flexible — any hierarchical pattern of fermion masses can 
be described by an appropriate AFS. 

The result (!24l) . which involves no flavor symmetry, has some similarities with the form (!25|) . 
First notice that a^, 6j, and Ck can be both positive and negative and therefore the factor of aibj 
(and Ckhj) in ^IM is statistically neutral. Anticipating this statistical averaging, we cast 
into the form of (1251) with the identification 
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e9 = e-^, e^ = e"5^. (26) 



An important feature is that the AFS factor is shared by all elements of both the up-type and 
the down-type Yukawa couplings that involve the left-handed quark doublet qj. This introduces 
a correlation between the up-type and down-type Yukawa matrices. 

We first study the probability distribution for a single entry in the Yukawa matrix, ignoring 
correlations with other entries. This allows us to determine the analogue of (jl]) for this Gaussian 
landscape. Here we do not have a distinction between the up-type and down-type Yukawa 
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Figure 4: Distributions of 10^ Yukawa matrix elements, generated from the Gaussian landscape 
on 5*^. From left to right the panels correspond to {d/L,gcs) = (0.08, 1), (0.1, 1), and (0.14, 1). 



matrices because we assume that the center coordinates of both the u and d wavefunctions 
are distributed randomly over S^. The up-type Yukawa coupling in ( 12^ is a function on the 
two-dimensional parameter space (a, 6) G [— L/2,L/2] x [— L/2,L/2]. The probability that the 
Yukawa coupling is larger than some value Aq is proportional to the area of an ellipse. 



^ 3(a-6)2 
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(27) 



When the signs of a and b are opposite, the overlap of the three wavefunctions is small and the 
Yukawa coupling becomes small. This is why the region of A > Aq is short in the (a — b) axis and 
long in the (a + b) axis. Comparing the area (1271) to the total area L^, we see the probability 
that A > Ao is given by 



P(A > Ao 

and hence the distribution is flat: 

dPjX) 
d\nX 
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(29) 



The distribution may cease to be fiat as a or 6 approaches ±L/2, where the parameter space 
ends, because we ignored the periodic boundary condition in the calculation that led to (!24|) . 
Setting this point aside we see that the probability distribution of Yukawa couplings in this toy 
landscape is fiat on the In A axis and has an approximate span of 



A In A = In 



A. 
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(30) 



arising from the inverse of the height of the distribution function (1291) . The overall hierarchy 
among Yukawa couplings AlnA = ln(Amax/Amin) is proportional to (L/d)^; the narrower the 
wavefunctions become, the further the wavefunctions can be separated, and the smaller the 
Yukawa couplings can be. As seen from (|2^ . the upper end of this scale- invariant distribution 
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Iogio(2^i2/7r) 



logio(2^^23/vr) 



logio sin 6^13 



Figure 5: Distributions of three quark Yukawa eigenvalues and mixing angles, based on a nu- 
merical simulation of the Gaussian landscape on S^, using {d/L,gcs) = (0.08,0.2). 



Amax is roughly ges\/L/d. Note that the scale- invariant distribution (j4]) was introduced almost 
by hand in [TU] in order to account for the large hierarchy among Yukawa couplings. It is 
interesting that this distribution is a natural prediction of our simple Gaussian landscape. 

We performed a numerical calculation to confirm the semi-analytical analysis above, taking 
account of the compactness of by making the wavefunction fl20|) periodic. The center coordi- 
nates of the quark wavefunctions ai,bj,Ck were generated randomly 10^ times, and the Yukawa 
coupling was calculated through fl2^ . This process was repeated for three different sets of the 
{d/L,gcs) parameters: (0.08,1), (0.10,1) and (0.14,1). The resulting distributions, shown in 
Figure m are all roughly scale-invariant (flat on a logarithmic scale), with heights proportional 
to {d/L)^, just as we expected from the semi-analytical discussion. 

3.2 Quark-Sector Phenomenology of the Gaussian Landscape 

Let us now study the probability distributions of the mass eigenvalues and mixing angles. Fig- 
ure O shows the result of a numerical simulation with {d/L,gcs) = (0.08,0.2). The distributions 
of Yukawa eigenvalues in Figure O are similar to those in Figure [H but with narrower distri- 
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10 si'^ ^12,23,13 

Figure 6: Distributions of the three AFS suppression factors, the three Yukawa eigenvalues, and 
the CKM mixing angles; the latter calculated naively from (15^ and (I5B1 - HU]) using the Gaussian 
landscape on S^. The CKM mixing angle 623 is the one most sharply peaked at zero, while ^13 
is most spread out. These figures correspond to d/L = 0.08 and thus Alog^ge = 5.65. Recall 
that these distribution functions are not reliable for small eigenvalues. 



butions of A^^^ and with the distributions of Ac,s shifted downward. The prominent difference 
between the mixing-angle distributions in Figure [5] and those in Figure [3] is the absence of the 
unwanted peaks at % ~ 7i/2. Thus we find the generation structure of the quark sector fol- 
lows from this Gaussian landscape; introducing correlations between the up-type and down-type 
Yukawa matrix elements works perfectly. Moreover, in contrast to Figure E] the distribution of 
6'i3 in Figure E] has a clear peak at ^13 <^ when displayed on a logarithmic scale. 

The distributions of Yukawa eigenvalues and mixing angles in Figure [5] can be understood 
analytically if we allow ourselves to make an approximation. We have seen in section 13.11 that 
both the up-type and down-type Yukawa matrices have an AFS structure. Thus, we begin by 
determining the probability distribution of the AFS suppression factors e"^'"'*^ in (l26l) . The value 
of e'^ (^g"''^) is determined by the distance |6| (|a|, |c|) of the left-handed (right-handed) quark 
wavefunction from the Higgs boson wavefunction. Since the center coordinates are scanned 
randomly over the extra dimension S^, the probability measure is 

dP{b) = yd\b\ , for < |6| < ^. (31) 

The measure for the right-handed quarks is the same, and we only deal with the left-handed 
quarks hereafter. Converting the variable from |6| to Ine"^ using (l26i) . we find 

dPiy) = 7rF, for < y < 1 , (32) 

where y = lne'?/Alne and we have defined Alne = -(L/2d)V3 = -{L/df/l2. 

The center coordinates of the three left-handed quark wavefunctions are chosen randomly, 
thus three AFS suppression factors follow (1321) independently. The smallest of these corresponds 
to the suppression factor ei for the lightest quark, while the middle factor €2 and largest factor 
€3 correspond to the middle and the heaviest quarks. The distribution of 61,2,3 is given by 

dP{yi, y2, 2/3) = j_ ^y"^ Q(^y^ _ ^2) e(t/2 - 1/3) , for < < 1 , (33) 
2^ ^ymys 
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where Ui = Inej/Alne. Note that yi > y2 > ys- The distributions of the individual AFS 
suppression factors €1^2,3 are obtained by integrating (15^ with respect to the other two variables: 

dP{y3) = l^-^—^^dy,, dP{y2)=3{l-^2)dy2, dP{y,) = l^.dy, . (34) 



These distribution functions are shown in Figure [6l Meanwhile, their averages are given by 

^^^^3) , \ (Inea) , . (Inei) (ln(ei/e2)) . - 

-r-j — = ?/3 =0.1, -r-j — = ?/2 =0.3, — — = =0.6, ———— = 1.5. 35 
A me Alne A me (In(e2/e3)) 

In a sense, this Gaussian landscape predicts the ratio of the AFS charges for the three generations: 
6:3:1. However this ratio only describes the statistical ensemble, and the distribution functions 
(IMl) contain more information. 

The distributions of Yukawa eigenvalues follow from (l34l) with the approximation 

ln(A,/A^ax)~lneM. (36) 

Explicit expressions of the distribution functions derived in this way are found in the appendix, 
fll84H186l) , and are plotted in Figure [61 From fll84H186l) we see that for small values of Zi = 
ln(Ai/Amax)/Alne (these correspond to large eigenvalues), the distribution functions behave as 
7^ 0, oc Z2, and oc zf, which is confirmed in the numerical simulation in Figure O Note that 
the distribution of A^^c begins at around 10"^, and that of X^^d at around 10~^, contrary to the 
approximate analytic distribution functions fll85lll86p . which begin at 0{1). This is probably 
due to the effects of diagonalizing the Yukawa matrices. As is familiar in quantum mechanics, 
two degenerate eigenvalues split even in the presence of the slightest perturbation. We later 
refer to this effect as the "diagonalization effect." 

We should also note that the AFS approximation (12^ [26]) breaks down for small values of 
Yukawa eigenvalues. The distributions based on the AFS approximation extend all the way 
down to ln(A/Amax) ~ 2Alne, while the numerical results cover a logarithmic range closer to 
fl30p . Since Alne ~ —A In A, the analytically derived range of In A covers twice the logarithmic 
scale that we expect. This discrepancy arises from the compactness of the extra dimension. 
That is, what really matters in the exponent of (124]) is 

mm [{a + + (b + ttiLY — {a + nL){b + mL)~\ , for n,mGZ. (37) 

As the center coordinates a and b approach ±L/2, non-zero choices of n and m may become 
just as important as n = m = in fl2^ . Indeed, when \a\ ~ \b\ ~ L/2 integers n and m 
can be chosen so that the last term is negative. When the compactness of the extra dimension 
is taken into account, the full expression 0371) cannot be larger than (L/2)^. This is why the 
distributions of Yukawa couplings and eigenvalues in Figure S] and E] span over ln(A/Amax) ~ 
Alne in the numerical results. The AFS approximation breaks down because the coefficient 
gij = e'^"'^^^^^'^^^^'^^'^/^'^^ is not statistically neutral when \ai\ and \bj\ are around L/2. 
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Distribution functions of the mixing angles can also be obtained by pursuing the AFS ap- 



proximation, along with the additional (crude) approximations: 

InVus ~ In (max{(e^/e^) 

u-sector; l/^2)d-scctor }) (38) 

\nVcb ~ In (max{(e^/e^) 

u- sector ; (el/4)d-scctoJ) (39) 
In Kb ~ In (max{(eVe3)u-sector, (e^ADd-sector}) • (40) 



Explicit expressions are found in fll88H190"l) and are plotted in Figure [61 These approximate 
analytic distribution functions capture qualitative features of the numerical results. Note that 
the distribution function of sin 6*13 becomes zero at logj^o sin 6*13 ~ because sin 613 can be of order 
unity only when all three eigenvalues are almost degenerate in either the up- or down-sector, 
and the probability for this to occur is small. The averages of the mixing angles in a logarithm 
scale are ordered 

(^13) < (^^12) < (^23) , (41) 

both in the numerical simulation and in the analytic distributions; see fll9ip . This is regarded 
as a consequence of the AFS charges in fl35|) : indeed in fl35|) we have 

(ln(6i/e3)) < (ln(ei/e2)) < (In(e2/e3)) . (42) 

Whether the assignments in (135|) are observationally acceptable or not is debatable, and we will 
return to this issue in section |H We discuss in section [5] how the "charge assignment" changes 
when the geometry 5*^ is replaced by other geometries. 

To summarize, we see that the qualitative expectations from an AFS-type mass matrix hold 
true in this landscape. In particular, the similarities between the landscape generated Yukawa 
couplings (124]) and those of (!25|) allow us to understand Gaussian landscapes, at least to some 
degree, using intuition based on models of AFS. Of course, between these approaches the origin 
of small parameters is completely different: in the landscape they arise from small overlaps 
of wavefunctions well-separated in the extra dimension, whereas in AFS they arise from small 
symmetry breaking parameters. Let us now emphasize this distinction. 

A crucial general feature of all AFS models is that the mass hierarchy between generations, 
^ m2 ^ nil, arises because there is a hierarchy in the amount of symmetry breaking coupled 
to these generations. This is true in the general Abelian case by the choice e-s ^ 62 ^ ei. 
In more restricted versions having ~ e*^% the hierarchy is imposed by a choice of charges 
Qi > Q2 > Qs- If the flavor symmetry is non- Abelian, then there is a hierarchy of symmetry 
breaking, for example the rank may be broken from z to z — 1 with strength e^. Thus AFS 
can in principle describe any flavor pattern, for example one heavy generation with Yukawas 
of order unity and two very light generations with Yukawas of order 10~^°. The situation with 
this Gaussian landscape is very different. Each Yukawa coupling is a statistical quantity, with a 
probability distribution that is approximately scale invariant over a range determined by a single 
small parameter, d/L. Relative to this range, the hierarchy of mass eigenvalues, including the 
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typical inter-generational mass ratios and mixing angles, arises purely from statistics. Unlike 
with the AFS parameters e^, there is no sense in which the fundamental theory distinguishes 
among generations. Therefore unlike with AFS, this Gaussian landscape cannot accommodate 
one heavy generation and two very light generations of comparable mass. Within statistical 
uncertainties, the landscape determines the AFS charges. 

3.3 Environmental Selection Effects 

It is a formidable task to understand all of the environmental effects that would be associated 
with a landscape scanning over the flavor parameters of the Standard Model. Furthermore, 
without a specific theory for the landscape it is unclear whether certain qualitative features of 
the flavor sector arise due to environmental selection, due to systematic features of the landscape 
distributions, or due to accident. Consider these examples. In the Gaussian landscape of this 
section, A1/A2 tends to be smaller than A2/A3. Therefore if d/L is chosen so as to explain 
the hierarchy A2/A3, then the relative lightness of the up and down quarks is explained. The 
value oi d/L required to reproduce the quark masses we measure may be selected dynamically 
within the fundamental theory or it may be selected due to environmental effects associated with 
having light up and down quarks. Likewise, the unexpected hierarchy mt/rrih may be due to 
the dynamical or accidental selection of different effective coupling constant (Jeff's for the up and 
down sectorsl^ or environmental selection could favor a very large top mass as described in [8]. 

We now describe qualitative effects associated with one possibility for environmental selection, 
which is the selection of a large top mass to ensure the stability of our present Higgs phase [8]. 
Consider a simple cut on the top Yukawa: 

logio \>- 0.3 . (43) 

We emphasize that we impose this t-cut to study qualitative effects; it is not intended to be 
precise. We first study the distributions of the quark Yukawa eigenvalues, which follow the 
distributions shown in Figure [3 The three distributions in the first row correspond to Yukawa 
eigenvalues before the environmental cut is imposed. Since our landscape has not introduced 
any difference between the up- and down-type sectors, the three distributions are the same in 
both sectors. Imposing the t-cut f l43l) reduces the sample to 5% of its original size. After the 
t-cut is imposed, the eigenvalue distributions are modified into those displayed in the second 
and third rows of Figure [71 A notable effect of the t-cut is that the distributions of the other 
Yukawa eigenvalues are dragged upward. This effect is more evident in the up-sector than in 
the down-sector, improving the fit of Xb/Xt, Xc/Xb, and Xs/Xc- 

One might consider that this improvement is not enough. Certainly the Xb/Xt prediction is 
improved, but we find that the Standard Model value only moves from the 3rd to about the 

^This hierarchy may also be due to weak-scale supersymmetry with a large tan/?; however this requires that 
the wavefunctions for the up-type and down-type Higgs are located at the same position in the extra dimension. 
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Figure 7: Distributions of quark Yukawa eigenvalues based on the Gaussian landscape on 5*^. 
The first row shows the distributions from Figure [5] before the t-cut is imposed. The last two 
rows show the distributions of the roughly 5% of Yukawa matrices that pass the t-cut. 

6th percentile. Furthermore, the typical value of A;, is large even without the t-cut, and the 
distribution of A;, is shifted upward after this cut is imposed. There are a number of ways to 
modify the Gaussian landscape of this section to alleviate this problem. Some of these have 
already been mentioned; others are presented in section \5\, where we study how the geometry of 
the extra dimensions affects the distributions of observables in Gaussian landscapes. 

The effects of the t-cut are also seen in the distributions of the mixing angles; see Figure [HI 
The probability that 623 ~ C>{1) is reduced significantly, which is understandable because A* (and 
hence A;,) tends to be very close to 0{1) in the sample passing the t-cut, while Ac and A^ are 
rarely larger than 10"^. The 6*23 distribution is reduced at small values as well, probably because 
the distributions of Ac and A^ are reduced at small values. Thus environmental selection on some 
flavor parameters modifies the distributions of other flavor parameters through a complicated 
chain of correlations. 
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Figure 8: Distributions of CKM mixing angles, based on the Gaussian landscape on S^. The 
first row shows the distributions before the t-cut is imposed, while the second row shows the 
distributions corresponding to the roughly 5% of Yukawa matrices that pass the t-cut. 

3.4 Summary 

We introduce a simple Gaussian landscape, based on a more microscopic description of the 
origin of fiavor involving a single extra dimension, which predicts probability distributions for 
the CP-conserving quark sector fiavor parameters. Using only two free parameters, ^^eff and d/L, 
this restrictive theory provides as good a fit to the fiavor parameters as conventional theories 
of flavor based on fiavor symmetries. This is especially so when the Gaussian landscape is 
augmented by environmental selection in favor of large top Yukawa coupling. This Gaussian 
landscape is characterized by the homogeneous scanning of the center coordinates of Gaussian 
wavefunctions in the extra dimension, which results in homogeneous scanning of Yukawa matrix 
elements on a logarithmic scale. The three eigenvalues of each Yukawa matrix then tend to be 
well separated from one another on the logarithmic scale. Thus, the hierarchical structure of 
quark mass eigenvalues is obtained. The wavefunctions of the three quark doublets qj and the 
Higgs boson h are assumed to be quite localized in the extra dimension, and the overlap of these 
wavefunctions introduces correlation between the up-type and down-type Yukawa matrices. This 
is the essence behind (the microscopic explanation) for the generation structure of the quark 
sector, that is the phenomenon that the heaviest (middle and lightest) up-type quark is coupled 
to the heaviest (middle and lightest) down-type quark in the 14^-boson current. 
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4 Goodness-of-fit Tests of Landscape Models 



Before proceeding to generalize the Gaussian landscape, we here take an aside to discuss some 
possibilities for testing theories that assume a huge landscape of vacua. A landscape theory that 
contains the observed vacuum is not directly falsifiable in the convention sense, since we can only 
compare one value of any "variable" (the value in our universe) to its predicted distribution. On 
the other hand, with a statistical theoretical prediction we can at least calculate the probability 
to have measured a more atypical value for a landscape variable. If this probability is small, 
the landscape can be excluded at a certain confidence level. Indeed, the more numerous are the 
predictions of a landscape theory, the greater the possibility to exclude the theory. Consider for 
example quantum theory, where statistical theoretical predictions are well-tested using scattering 
experiments that count the event rates in a large number of bins/modes. The parameterization 
of the Standard Model is dominated by the 22 flavor parameters, 20 of which have either been 
measured or may be measured in the near future. Therefore it seems that the flavor structure 
of the Standard Model provides the best realm in which to test landscape theories. 

In this section we study two statistical tests, the chi-square statistic and the p- value statistic, 
to illustrate how Gaussian landscapes may be excluded by experiment. As this section is intended 
to be only illustrative, we consider only the quark sector as described by the Gaussian landscape 
on (which in this section is referred to as the model) and we do not consider any weight 
factors that may arise due to cosmological evolution or environmental selection effects. The 
existence of correlations among the distributions of flavor parameters in Gaussian landscapes 
allows for a rich statistical analysis, and the discussion below only scratches the surface. 

4.1 The Chi-Square Statistic 

Our first example is to use the so-called chi-square statistic as a measure of goodness-of-fit 
between the observed flavor parameters and their hypothetical landscape distributions. The 
CP-conserving flavor parameters in the quark sector consist of the nine variables: 

{Xj} = logio{A„, Ac, \t, Ad, A^, Afe, 6^12, 623, sin6'i3}, for J = 1, ■ ■ ■ ,9, (44) 

with the logarithm acting on every element of the list. The CKM phase is not included in this 
list because the model we consider here is CP-conserving. The chi-square statistic is 

= E (^^ - ) i^'')iJ (^^ - ) ' (45) 
I, J 

where hats denote the measured value of a parameter. The values listed in Table [T] are used for 
{Xj}. Brackets denote the landscape average, and the co-variance matrix V is given by 

Vij = Cov{Xj, Xj) = {XjXj) - {Xj) {Xj) . (46) 

The chi-square statistic is invariant under any linear transformation among the variables {Xj}. 
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The correlation among the parameters in fH4j) cannot be ignored. Specifically, for the three 
mass eigenvalues in the up sector, the correlation matrix, p/j = Vjj I \JViiVjj, is given by 

/ 1. 0.57 0.39 \ 
Pjj=\ 0.57 1. 0.58 , for /, J = 1,2,3. (47) 
\ 0.39 0.58 1. / 

Note that the off-diagonal terms are not negligible compared to the diagonal terms. The pos- 
itive correlations are as expected from the combined probability distribution (133|) of the AFS 
suppression factors log^Q ef Meanwhile, the correlation between the up and down sector mass 
eigenvalues, {log^o A„, log^o A^, log^o A*} x {log^o A^, log^o A,, log^o AJ, is given by 

/ 0.32 0.27 0.20 \ 
Pij=\ 0.27 0.39 0.29 , for J = 1, 2, 3 and J = 4, 5, 6 . (48) 
\ 0.20 0.29 0.43 / 

This is also sizable and positive, confirming that the mass eigenvalues in the down sector are 
dragged upward under the t-cut. The mixing angles are also correlated with mass hierarchy; the 
correlation matrix between log]^o{A«/Ac, Ac/At, Au/A^} and log]^o{^i2) ^23, sin 6'i3} is given by 

/ 0.34 -0.14 0.28 \ 
pjj = -0.21 0.42 0.14 , (49) 
\ 0.16 0.21 0.38 / 

where now I = Xi — -^2 ~ -^1 ~ -^3 and J = 7, 8, 9. We see that larger mass hierarchies 
are correlated with smaller mixing angles. 

The covariance matrix Vu determines the principal axes Yi = CijXj, with the covariance 
matrix Vij being diagonal when the basis {Xj} is changed to {Yi}. Since Vij is diagonal the 
observables {Yj} are independent at least up to second order. If the distributions of the li's 
were completely independent, and Gaussian, then the quantity would follow the chi-square 
distribution for random measurements of [Yi = CijXj}. This distribution has probability density 

4^ = 2^/2r(iV/2) x^'-'e-^"^' , (50) 
dx 

where is the number of independent Gaussian random variables entering (H5|) . This distribu- 
tion has a mean of and a standard deviation of \/2N. Thus the prediction x^ ~ A^(l ± a/2/A^) 
is very sharp when a large number of observables are predicted by the landscape. Note that too 
large a value of indicates that the observed fiavor parameters are atypical of what is expected 
from the landscape distributions, while too small a value of indicates that the correct un- 
derlying theory has less randomness than is exhibited by the landscape in question. Therefore, 
even though we can measure only one value for each observable, and hence measure only one 
value of x^, the chi-square statistic can still be a powerful tool in testing landscape theories. 
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dP/dx^ HdP/dx^) 

Figure 9: The left panel displays the probability distribution of obtained numerically (his- 
togram) using data generated by the model, along with the actual chi-square distribution for 
N = 9. The right panel displays the logarithm of the numerical distribution of (histogram) 
and a normal distribution with a mean = 9 and standard deviation \^2N = \/T8. 



Our model predicts distributions for nine flavor parameters in the quark sector, but it 
contains np = 2 free parameters (gcs and d/L) that can be fixed by hand and tuneco to fit 
two out of nine of the l^'s, rendering two terms in Ylii^i ~ vanish. Thus we have 

N = 9 — Up = 7 independent degrees of freedom. Therefore if the distributions describing the 
Yi were Gaussian, we would expect x^ ^ 7 ± 4. In fact the model does not predict Gaussian 
distributions for any of the flavor parameters. However, Figure IH] shows that the distribution of 
X^ f l45l) using data generated by the model looks quite similar to the corresponding chi-square 
distribution, which has N = 9. Furthermore, the central limit theorem guarantees that as is 
increased the distribution of x^ approaches a normal distribution with a mean N and a variance 
that grows in proportion to A^. can be as large as 20 for the observable flavor parameters of 
the Standard Model, assuming a landscape with at least 20 scanning parameters. 

We calculate the chi-square statistic for the nine CP-conserving quark-sector flavor parame- 
ters using the 5*^ model to calculate (Xj) and V/j, and using Tabled] for Xj. This gives x^ ~ 6, 
which agrees very well with the prediction x^ ~ 7 ± 4. That the obtained value is neither too 
large nor too small indicates that the observed flavor parameters are typical and that they ex- 
hibit the expected level of statistical fluctuation in this landscape. The same analysis can be 
carried out for the 5*^ model after a hypothetical environmental weight factor favoring a large 
top mass has been imposed. Using only model data that pass the t-cut (l43l) . we find x^ ~ 8, 
which also agrees well with the prediction. Imposing the t-cut increases the value of x^? but 
one should keep in mind that the parameters g^s and d/L were chosen to fit the observed flavor 
parameters before imposing the t-cut, without respect for the ensemble that passes the t-cut. 

"'^"In fact the free parameters g^s and d/L have not been fully exploited, since we have not performed a maximum 
likelihood analysis to determine the values that minimize in ()45p . However, the parameter values we use have 
been chosen to qualitatively fit our expectations for the measured values of {^/}- 
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4.2 The P-Value Statistic 



The chi-square statistic has a disadvantage as well. Although the central limit theorem guar- 
antees that for large that the distribution of ^ approaches a normal distribution, regardless 
of most of the details of the actual distributions of the {X/}, is not more than 20 for the 
application to flavor parameters. Moreover, some of distributions of the flavor observables have 
sharp cut-offs, such that the probability to obtain a large value of [Xi — {Xi)Y may be much 
smaller than what is expected from a normal distribution with a mean (X/) and a variance V//. 
The chi-square statistic does not completely account for such a situation. 

On the other hand, the p-value statistic is capable of handling non-Gaussian distribution 
functions. The p-value of a single flavor parameter Xj is simply the fraction of the distribution 
that is more atypical than the measured value X/. Specifically, we may consider the p-value 
to be the fraction of the distribution that is farther from the median Xi than is X/. If we 
approximate the probabihty density as continuous with profile f{Xi), then the p-value pi is 

/X[ /-oo 
2fiXj)dXj hi Xi<Xi, pi{Xi)= 2f{Xi)dXi forX, >X,. (51) 
oo JXi 

Note that like x^? Pi is a random variable for every independent measurement X/. However, 
whereas is distributed according to the chi-square distribution, pi is distributed uniformly 
between zero and one. The p-value is also the probability that the Kolmogorov-Smirnov statistic 
(the value of the D function) of one sampling of the variable X/ could have been larger than 
what is calculated from X/. Hence a measurement giving p-value pi excludes the corresponding 
theory at a confidence level of max{p/, 1 — pi} x 100%O 

Reference [15] used the p-value statistic to test the hypothesis of Neutrino Anarchy. There, 
the probability distribution of the lepton sector flavor parameters is factorized into those of 
mass eigenvalues, of each individual mixing angle, and of the CP phase [13]. Then the p-value 
statistic can be applied separately to each of the three mixing angles. In the case where the 
statistical variables X/ are independent of each other, k = YliPii-^i) represents the fraction of 
the distribution of {X/} for which each element Xj is more atypical than the measured value 
X/. The probability that k = YliPii-^i) could have been smaller than k is given by (e.g. [T5] ) 



P 



n^^^=^E^(-i^^) • (52) 

/ a=0 ■ 



Here a is simply a summation index used to simplify the last expression, and N is the number 
of degrees of freedom in {X/}. For a random measurement of a set of independent variables 
{X/}, the p-value calculated through (!52|) is distributed uniformly from zero to one. Thus a 



^^P-valucs very near zero indicate that the measured values {Xj} are collectively atypical of their predicted 
distributions; p- values very near unity indicate that the measured values do not exhibit the randomness expected 
from the predicted distributions. 
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random measurement giving p-value p can exclude a landscape theory with a confidence level 
max{p, 1 — p} X 100%. If we treat the nine CP-conserving fiavor observables in the quark sector 
as if they were independent, then the data Xj and the hypothetical distributions dP[Xi) / dXj 
of the model can be used to calculate the p-value through Here, = 7 is used in ^% 
because the values of and d/ L can be tuned to bring two of the X/ to be the same as their 
medians. We find p ~ 0.66, suggesting a good fit with observation. 

In fact, however, the statistical variables {Xj} are not independent. To the best of our 
knowledge, there is no standard definition of the p-value statistic for correlated multi-variable 
distributions. One possibility to account for correlation is to calculate k by actually counting the 
fraction of randomly generated parameter sets {X/} that have each Xi more atypical than that 
of the measured set {Xi}. This is as opposed to using k = Y[jPi{Xi), which gives this fraction 
in the absence of correlations. In a sample of 5 x 10^ data sets generated using the numerical 
simulation of the model, we find 52 such sets, giving k = 2^x 

(5^) and hence p ~ 0.73±0.03. 
The uncertainty is obtained by assuming the variance in counting A'^atyp atypical data points is 
of the order of A'^atyp- Yet another attempt at a "p-value like" measure of the goodness-of-fit 
may be to judge atypicality by using in (15T|) the distributions of the principal axes Yi, their 
measured values Yi, and their medians Yi = CijXj — as opposed to using X/, X/, and Xj — since 
the statistical variables Yi do not have correlations, at least up to second order. The fraction of 
{Yi} more atypical than the measured set {Yi} can be calculated using k = Y[iPi(Xi)'^ using ( l52ll 
we find this gives p ^ 0.42. If there were truly no correlations among the Yi this would give the 
same result as calculating k by actually counting the fraction of randomly generated parameter 
sets {Yi} that have each Yi more atypical than that of the measured set {Yi}. However, in our 
ensemble of 5 x 10^ data sets we find just one for which each element is more atypical than its 
corresponding measured value, giving p ^ 0.19;']o;59. 



5 Geometry Dependence 

In section [3] we introduced and analyzed a Gaussian landscape based on a single extra dimension, 
and found that it could provide the hierarchy, pairing, and generation structures of the observed 
quark sector. However, the number and geometry of extra dimensions need not correspond 
to the Gaussian landscape. Therefore, we initiate a study into how the compactification 
geometry of extra dimensions affects the probability distributions of observables. We find that 
the qualitative results of section [3] can be achieved by Gaussian landscapes in other geometries of 
extra dimensions and that the distribution functions of fiavor observables are largely insensitive 
to the details of these geometries. However, the behavior of the distribution functions of mass 
eigenvalues at large values is affected by the number of extra dimensions. In section [TJ we argue 
that the geometries used in Gaussian landscapes are related to base manifolds of torus-fibered 
geometries used in the compactification of string theory. 
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Figure 10: Distributions of 10^ Yukawa matrix elements, generated from the Gaussian landscape 
on T^. From left to right the panels correspond to {d/L,ges) = (0.1, 1), (0.12, 1) and (0.14, 1). 

5.1 Gaussian Landscapes on and S'^ 

To investigate the robustness of the Gaussian landscape on S^, we study other geometries. We 
first look at the simplest extension to more than one dimension, T"^ = x S^, and focus on 
a "square torus" where the two periods of the torus are both L, and the directions of the two 
are orthogonal. We assume that each of the quarks and the Higgs have rotation-symmetric 
Gaussian wavefunctions of the form 

I y~yo I ^ 

^{y;yo)o(.e (53) 

where the center coordinates of each particle are randomly scanned over the internal space 
T^. The up-type and down-type Yukawa matrices are calculated by the overlap integration fl22|) . 
which is naturally generalized to integration on T^. We defer to section [7| a discussion of the 
extent to which these assumptions result from a dynamical field theory on extra dimensional 
spacetime, and for the moment focus on the phenomenology of this landscape. 

Figure [TU] shows the distribution of Yukawa matrix elements for a numerical simulation 
of this model. As in the model, the only parameters relevant to these distributions are 
9cS = 9/{MqL) and d/L, where now Mq is the cut-off scale of the effective theory in 5+1 
dimensions. Furthermore, a larger hierarchy is generated when the wavefunctions are more 
localized, i.e. when d/L is smaller. The key difference from is that the distribution of Yukawa 
matrix elements is no longer scale invariant. Instead, on the probability density for the largest 
and the smallest matrix elements is depleted. With regard to the Yukawa eigenvalues and mixing 
angles, we find that a hierarchical pattern of Yukawa eigenvalues is generated (Figure fTT]) and that 
the distributions of mixing angles are peaked at 9ij = but not at 9^ = 7r/2 (Figure [T^. Thus 
the flavor structure of the quark sector follows from the Gaussian landscape on T^. Despite the 
apparent difference between the distribution of the Yukawa matrix elements of the two landscapes 
(comparing Figure H] and Figure [TO!) . we see that the distribution of masses and mixing angles are 
roughly the same when we compare Figures [11] and [12] to Figures [7] and [8] In particular, these 
distributions all come with a width of about an order of magnitude, and the differences between 
the distributions from the two toy landscapes is not, statistically, very significant compared to 
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Figure 11: Distributions of Yukawa eigenvalues in the Gaussian landscape on T^, based on a 
numerical simulation with {d/L,gcs) = (0.1,0.4). The first row shows the distribution of the 
three eigenvalues of the up (and down) sector. The second and third rows display the eigenvalues 
of the roughly 6% of matrices that survive the t-cut of section [231 



this width. This demonstrates that the flavor structure in the quark sector is a robust feature 
of Gaussian landscapes. 

The biggest difference between the distribution of observables in the two toy landscapes is 
a more left/right symmetric probability distribution for logio A^^f logio As, c on T^. This is a 
consequence of the difference in the distribution of individual Yukawa matrix elements, Figure [TU] 
vs. Figure |H The depleted probability density of the largest Yukawa matrix elements in the 
Gaussian landscape on results in reduced probability for largest values of A;,,* and Xs,c- The 
distribution of Xb^t becomes much broader and less peaked at the largest possible value, Amax- 
This allows the bottom Yukawa coupling to become as small as its measured value rather easily, 
while still accommodating the measured value of the top Yukawa coupling. 

The effects of a possible environmental selection for a large top Yukawa coupling are studied 
in Figure [TT] and in Figure [T5], using the cut condition as a crude approximation to the 
environmental selection effect. The distributions of the mixing angles 623 and ^13 are shifted 
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Figure 12: Distributions of CKM mixing angles in the Gaussian landscape on T^, based on a 
numerical simulation with {d/L, Qcs) = (0.1, 0.4). The bottom row displays the mixing angles of 
the roughly 6% of matrices that survive the t-cut of section 13.31 



toward smaller angles, just as on S^. Since the cut condition is in favor of a larger hierarchy 
between the lighter quarks and the heaviest quarks, smaller 623 and 6*13 are natural consequences. 
The inequality (Ac) > (A^) also follows from the cut, just like in section [31 The most impor- 
tant difference between the two landscapes may be in the distribution of the bottom Yukawa 
coupling after the cut is imposed. It was rare that A;, be less than 10~^ in the lower-right 
distribution of Figure [71 but a significant fraction is below 10^^ in Figure [111 Therefore the 
observed hierarchy Aj/Af, may be understood within the context of a Gaussian landscape with 
an appropriately chosen geometry, especially when there is environmental selection for a large 
top Yukawa coupling. 

For comparison, we also perform a numerical study of a Gaussian landscape defined on 
S"^. As before, the quarks and Higgs are represented by localized wavefunctions with width d. 
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Figure 13: Distributions of 10^ Yukawa matrix elements, generated from the Gaussian landscape 
on 5^. From left to right the panels correspond to {d/R,gcs) = (0.4, 1), (0.5, 1) and (0.6, 1). 



Specifically, on 5^ a wavefunction centered at ^ = is given b; 




V9(6',0) oc e ^wW . (54) 

Of course, to generate an ensemble of Yukawa matrices the central coordinates of each wave- 
function are scanned independently and uniformly over the geometry S*^. Note that like the 
previously described landscapes, the Gaussian landscape on S"^ is characterized by two free 
parameters, d/R and Qck = g / {^/tiMqR). 

Figure [13] shows the distribution of Yukawa matrix elements for this model. Note that the 
overall shapes of the distributions are remarkably similar to those in Figure [IHl which correspond 
to T^. Figure [H] displays the distributions of the three mass eigenvalues for the Gaussian 
landscapes on 5^ and T^. These distributions are quite similar between S"^ and T^. We find that 
the other phenomenological aspects of the Gaussian landscape, discussed above with respect to 
the geometry T^, hold true on S'^ as well. 

5.2 Dependence on Dimensionality and Geometry Independence 

Having seen that the hierarchical mass eigenvalues and the generation structure of the quark 
sector are robust predictions of Gaussian landscapes, the next subject of interest is to understand 
to what extent the details of distribution functions depend on the internal geometry. Already we 
have seen slight differences between the Gaussian landscape on 5"^ and the landscapes defined 
on and 5^, yet remarkable similarity between the distributions coming from and S*^. We 
would like to understand where these differences and similarities come from, without having to 
run simulations on all possible choices of internal geometry. To do this, we recall the analysis 
of sections 13.11 and 13. 2[ where the distributions of flavor observables on were understood 
analytically within the context of the AFS approximation. As we show in this section, this 
analysis can be generalized to Gaussian landscapes on any geometry of extra dimensionso 

-'^^ Although this wavefunction is not smooth at 6' = tt, this is not of present concern. That is, the purpose of 
this numerical simulation is not to determine the distribution precisely but to study its qualitative aspects. 

^^Tliere is also a practical motivation for this approach. The numerical integration time involved with per- 
forming overlap integrals on extra dimensions grows very large as the number of extra dimensions is increased; 



31 



-8 " 1 -2 -8 -6 ! _ -8 -6 - 



400p . . . . . . . ■ . . . . . . . ^ r ' ' ■ ' ' ' ■ ' ' ' ■ ' ' ' ■ n 600 




logio Kd logio loglO Kb 

Figure 14: Comparison between the distributions of Yukawa eigenvalues from the Gaussian 
landscape on with {d/L,gcs) = (0.1,0.4) (top row) and the Gaussian landscape on 5*^ with 
{d/R,g^s) = (0.45,0.4) (bottom row). 



Let us consider a D-dimensional internal space with local coordinates y. Then in the limit 
where we can ignore the finite size of the internal space the Yukawa coupling matrix is given by 

A^^. oc e-3^(l'^'l'+l^^^l'-'^'-^^\ (55) 

which replaces A*'^- in fl24l) for the Gaussian landscape on S^. Here Si and bj are the central 
coordinates of the quark wavefunctions, relative that of the Higgs. The Yukawa matrix again 
has the AFS form fl25l) . with suppression factors 

e'j = e "^2- ^ = e ^ , (56) 

replacing (126|) . In the extra factor gij = e°i-*j/3'^^ (which is not necessarily of order unity) the 
vectors and bj are sometimes unaligned, sometimes parallel, and sometimes anti-parallel. 
This generates a random coefficient to each Yukawa coupling that is statistically neutral in the 
AFS approximation. We note again that this analysis is valid only when the compactness of 
the internal space is unimportant and the local geometry can be approximated as a flat D- 
dimensional space. This is equivalent to focusing on only the largest Yukawa matrix elements. 
Ignoring the statistically neutral factor ^'^i-bj/^d'^ ^ Yukawa matrix elements are roughly 



r 



2 



A~e -lnA = ^, (57) 



meanwhile generating ensembles of flavor observables involves performing large numbers of these integrals. 
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with r = (a, b) scanning a 2D-dimensional space. The natural probabihty measure is 

dP^—j^ oc ' ' ' ' oc - (-lnA)^^^d|lnA|, (5^ 



where L is the typical size of the extra dimensions. Indeed, the distribution of the Gaussian 
landscape on S*^ (-D — 1 = 0) is flat, c.f. Figure HI and those of the Gaussian landscapes on 
and 5*^ (-D — 1 = 1) are linear in | In A| at their upper end0 These results are independent of the 
value of d/L. The logarithmic range of the distribution of Yukawa couplings scales as {L/dY- 

The distribution of the AFS suppression factors, e^, can be obtained as in section [31 In more 
than one dimension the measure (I3T!) is generalized to 

d^h . r.dV(b)db^ , , 

where h = \b\ and V{b) is the volume enclosed within a distance b from a given point. Using 
(|56|1 . it can be converted into a distribution of the AFS suppression factors. 



dPiy) = fiy)dy, 



Vtot db^ 



(60) 



Here, Ine is normalized by Alne = — |(&max/'^)^, so that y = Ine/Alne = (6/fcmax)^ ranges from 
zero to one. For example, 6max = L/2 for S^, bmax = for and 6max = ^ri? for S"^. The 

ratio (6max/ dY determines the overall logarithmic range of hierarchy, and the volume distribution 
function f{y) controls the shape of the distributions within y G [0, 1]. It is straightforward to 
find f{y) for a given geometry of extra dimensions; for example, for S^, S*^ and S^, 

f ( \ ^ f f \ TT sin(7r^) sin^(7ryy) 

/s^(y) = ^, fAy) = j^j=-, fs-4y) = ^^, (6i) 



and for the "square torus" T 



2 



fr^iy) = I - 2 arcsin 1 " ^ j " V2) • (62) 

These volume distribution functions are displayed in the first row of Figure [151 

The distribution of the smallest, middle, or largest AFS suppression factor is obtained by 
integrating the other two variables out of the probability distribution 

dPivi, 1/2, y-i) = 3!/(yi)/(?/2)/(y3)0(z/i - Z/2)6(?/2 - Z/s) dyidy2dyz , (63) 

where yi = Inej/Alne and we remind that in this notation yi > y2 > ys- The AFS sup- 
pression factors yi for S^'^'^ are displayed in the second row of Figure [13 The mass eigenval- 
ues A"*-'^'' are approximated by e1e^^'^\ and hence the distribution of Zi = ln(A"*^'^V'^max)/Alne 



^^Tlie distributions on are not as precisely linear for large matrix elements as those on T^. This can be 
understood in terms of geometry: 5*^ has positive curvature while is flat. 
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is given by a convolution of the distribution function of = Inej/Alne (third row of Fig- 
ure [T^ . The diagonahzation matrices for the up and down sectors are determined by three 
diagonahzation angles, which are approximately equal to ei/ej, i < j . Thus, the distribution of 
Aij = ln(ej/ej)/A Ine = Hi — yj is obtained by integrating the variables 1/1,2,3 out of (jHSD while 
keeping the distance A^j = yi — yj (fourth row of Figure [15]). Finally, the CKM mixing angles are 
obtained from the diagonahzation angles by approximating the mixing angles to be larger one 
of Aij in the up sector and the down sector. These are displayed in the fifth row of Figure [151 
Note that the analytical results for the mass eigenvalue distributions (the third row of Figure [T5l) 
capture the qualitative features of the numerical results for the Gaussian landscapes on S^, 
and S"^ very well. 

It is now possible to understand why the distributions of flavor observables are quite similar 
between the Gaussian landscapes on and S'^. All of these distributions are calculated using 
only the volume distribution function f{y). Although f{y) is quite different for these two geome- 
tries, as is seen in the first row of Figure [TSj a number of consecutive integrations is required to 
obtain distribution functions for observable flavor parameters. For any AFS suppression factor 
yi oc Inej, the volume distribution function f{yi) in (!63l) is not integrated, and thus singular- 
ities in the original f{yi) remain in the distribution of log^g^j- Nevertheless, the difference in 
f{y) between and S'^ is already less evident in the distributions of log^^Qe^, (second row of 
Figure [T^ . because the integration of two variables out of (1^ takes a certain average of f{y)- 
Meanwhile, to obtain the distribution functions of the mass eigenvalues, diagonahzation angles, 
and the CKM mixing angles, each f{yi) is integrated at least once, and the geometry-dependent 
information contained in f[y) is further smeared out. This is why the distributions of the Yukawa 
eigenvalues are quite similar between and S"^ in Figures [TH and [151 Since the distributions of 
CKM mixing angles involve a greater number of integrations, these distributions are even less 
sensitive to the underlying geometry. 

The analytical discussion so far explains why the distributions of flavor observables are similar 
for the Gaussian landscapes on and 5^, but it remains unseen why the distributions of mass 
eigenvalues of the Gaussian landscape are somewhat different from those of and S"^. 
Therefore, we now focus on how the number of dimensions in the internal geometry affects the 
distributions of flavor observables. As we have already seen in the discussion surrounding (!58|) . 
the number of extra dimensions directly affects the shape of the volume distribution function. 
Specifically, we find f(y) oc y^^'^~^ for y greater than but near zero. Thus for very small 

dP{yi) oc yl'^''^~^dyi , dP{y2) oc "M?/2 , dP{y^) oc y^^^'^dy^ . (64) 

Distribution functions of the masses and mixing angles can be obtained by using the approxi- 
mations (l36l) and (l38l - [40|) . analogous to the analysis earlier in this section. In the limit of large 
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Figure 15: The volume distribution functions f{y), AFS suppression factors logio 61,2,3, the 
Yukawa eigenvalues logio Ai, 2,3, diagonalization angles logio(ei/ej) {i < j), and the CKM 
mixing angles logio sin for S^, T^, and S^. In all cases the broadest mixing angle distribu- 
tion corresponds to ^13, while the distribution most sharply peaked at zero corresponds to 6*23. 
These use approximations ( |36l) and ( |38l - l40l) . which are not reliable for small Yukawa eigenval- 
ues. For clear comparison the variables yi, Zi, and ttj have been converted to observables using 
a common logarithmic scale, A logio e ~ —5.66, that corresponds to the scale of the Gaussian 
landscape on with d/L = 0.08. 



eigenvalues and large mixing angles these distribution functions behave as as 

(65) 



dPj^l) ^ ,3D-l dP{z2) 2D-1 dPjZs) n_i 

LX x^i , LX , ^ ^3 

azi az2 dzs 



^EM^o, 'JM^o, ^^oc*. 

dti2 dt23 dti3 



where again Zi = ln(Aj/Amax)/Alne and tij = In sin % /A In e. Note that the behavior of the 
distribution functions of the mixing angles does not depend on the number of dimensions. The in- 
dependence above is quahtatively confirmed in Figure [T3 Because of the slow rising of the = 
ln(At^b/Aniax)/Alne distribution function in higher dimensions, the weight of the distribution of 
^3 is shifted toward larger values, allowing for smaller values of A;,. 

As is evident from Figure [151 however, the approximation f{y) oc y^/^-i is valid for only 
a very narrow range of y near zero. Using the explicit form of f{y) instead, we can extract 
another systematic D-dependence. The statistical average of yi = Inej/Alne plays the role of 
an AFS charge of the ith generation, and yi — yj determines the hierarchy between the ith and 
jth generations. For sequences of geometries with different numbers of extra dimensions — S^, 
S"^, and 5*^, and = 5*^, T^, and = x x — an explicit calculation reveals 



This ratio clearly decreases in Gaussian landscapes with a greater number of extra dimensions, 
meaning that the middle eigenvalue becomes (statistically) closer to the smallest eigenvalue on a 
logarithmic scale. This also means that in in the AFS approximation the diagonalization angle 
(and hence the mixing angle) between the first and second generations becomes more significant 
compared with that between the second and third generations. This tendency is confirmed 
in Figure [15], where the 612 distribution becomes closer to the 623 distribution for larger D. 
Given that we measure Xd/^s > Xs/Xb and 612 > 623, this might be regarded as an indication 
that D > 1. Note however that the compactness of extra dimensions affects the distributions 
of the smallest eigenvalues, and therefore the distributions of observables based on the AFS 
approximation cannot reliably be used to infer the number of extra dimensions. 

Let us now summarize the conclusions of this section. The basic features of flavor in the 
Standard Model follow from Gaussian landscapes on any geometry of extra dimensions. Specif- 
ically, the distribution of Yukawa matrix elements is generally broad and thus there is always 
hierarchy among Yukawa mass eigenvalues. Furthermore, the up-type and down-type Yukawa 
matrices are correlated due to the common overlaps of left-handed quark doublets and the Higgs, 
so that generation structure is always realized. Even in the details of the distribution functions of 
flavor observables, there is not much geometry dependence. We understand this in terms of the 
AFS approximation, where all of the geometry dependence is encoded in the volume distribution 
function f{y). This function is integrated and convoluted many times to obtain the distributions 
of observables, thus smearing out the original geometry dependence. This means that we cannot 
learn very much about the geometry of extra dimensions from the observed masses and mixing 
angles. On the other hand, it appears that we can understand the qualitative pattern of masses 
and mixing angles without knowing much about the details of the underlying geometry. The 
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number of extra dimensions, however, leaves its footprint on some distribution functions near 
y ~ 0, because the boundary at y = remains a boundary in some convolutions. We also find 
that the mass eigenvalues of the second generation become closer (statistically) to those of the 
first generation as we increase the number of extra dimensions in Gaussian landscapes. 

For various geometries, the volume distribution functions f{y) are usually moderately varying 
functions of y over the range y G [0, 1]. Since the distribution functions of the AFS suppression 
factors yi oc Ine^, the mass eigenvalues Zi cx InAj, and the mixing angles tij oc In sin 6*^ are all 
derived from f{y), these distribution functions are moderately varying functions of the loga- 
rithmic variables yi, Zi, and tij. The widths of these statistical distributions are quite broad; to 
achieve otherwise would require an exponentially steep f{y), which does not happen in Gaussian 
landscapes. Therefore, the Gaussian landscapes in this article are different from those in |16j . 
where all of the dimensionless coupling constants of the Standard Model have narrow-width 
Gaussian distributions. This feature of our toy landscapes is traced back to our assumption that 
the localized wavefunctions become exponentially small as one moves away from the centers of 
localization. 

6 Including the Lepton Sector 

Yukawa couplings generated from simple Gaussian landscapes are in good qualitative agreement 
with the mass hierarchies and small mixing angles in the quark sector. On the other hand, the 
lepton sector is characterized by large mixing angles and very small neutrino masses. We now 
turn our attention to how these qualitatively distinct features might be explained within a single 
landscape model. With the challenges of describing the lepton sector comes the opportunity 
to predict the probability distributions of three yet-to-be measured observables in the lepton 
sector: the mixing angle 6'i3, the CP phase in neutrino oscillations and the mass parameter rnpis 
of neutrinoless double beta decay. 

In order to accommodate large mixing angles in the lepton sector, in section 16.11 the Gaus- 
sian landscape is extended to include different widths for wavefunctions of particles in different 
representations. We will see in section 16.41 that the mixing angles can be large indeed, and we 
also learn that complex (CP violating) phases play a crucial role in determining the distribution 
of mixing angles (and vice versa in section [6751) . Before studying mixing angles we study the 
charged lepton mass spectrum in section 16. 2[ Small neutrino masses are assumed to be gener- 
ated via the seesaw mechanism, and in section 16.31 we describe how to generate the statistics of 
the right-handed neutrino Majorana mass terms in Gaussian landscapes. Finally, in section 1^31 
we impose some cuts on the statistics generated by a particular version of the Gaussian land- 
scape. This is done to obtain a general feeling for how the probability distributions of yet-to-be 
measured observables can be affected when distributions conditional on measurements already 
performed in this universe are considered. 
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Note that our goal in this section is to find a single theoretical framework that can describe 
the various fiavor structures that are observed; i.e. we aim to identify what subset of landscapes 
share key qualitative features with the observed fiavor structure. Ultimately, the success of the 
landscape picture will depend on both the existence of phenomenologically viable subsets to the 
landscape and that these subsets are not too atypical of what is expected from the full landscape, 
after cosmological and environmental selection effects are accounted. 

6.1 Landscapes with Delocalized 5 = {d,l) 

In the traditional AFS approach, the large mixing angles of neutrino oscillation and a mild 
charged lepton mass hierarchy result if the three lepton doublets are not strongly distinguished 
by the fiavor symmetry. For example, an AFS may be broken by a single parameter, with 
the three g C 10 = (g, u, e)'s strongly distinguished by the symmetry charges while the three 
I G 5 = {d, /)'s are not. Since Gaussian landscapes can mimic an AFS structure in the Yukawa 
couplings, the idea of different strengths of flavor symmetry breaking for 5 and 10 can be 
translated into the framework of Gaussian landscapes. 

As is seen in sections E] and El the ratio of effective AFS charges for the three fermions in a 
given representation is determined (statistically) by the Gaussian landscape; we have no freedom 
to choose these by hand. On the other hand, the overall hierarchy depends on the parameter 
Alne oc — (L/rf)^. So far, for simplicity the parameter d has been chosen to be the same for 
all of the Standard Model wavef unctions; however this parameter can be different for flelds in 
different representations. If the 10 and 5 flelds have different width parameters d\o and d^, then 
the AFS suppression factors associated with these flelds are also different. In such a landscape 
the single parameter d is replaced by three: dn, dio and d^. 

If the wavefunctions of the fermions in the 5 representation are not particularly localized, the 
overlaps between the wavefunctions of the 5 flelds and the Higgs boson do not vary hierarchically 
with the peak locations of the 5 Gaussian wavefunctions. Therefore the Yukawa couplings 
associated with the three 5 flelds are not hierarchically separated when d^ is not much less than 
L. This is along the line of the idea in [17]. On the other hand, such a choice of maintains 
the hierarchical structure of the masses and mixing angles in the quark sector, when d\o/L and 
dn/L are chosen the same as before. The main exception is that the hierarchy of the down-type 
quark masses becomes a little smaller, as we see below, which is actually in good agreement with 
observation. 

6.2 Charged Lepton Mass Spectrum 

The masses of the charged leptons derive from the Yukawa interaction 

C = Xl,eakh\ (68) 
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where A^j is generated in analogy to fl22l) . except with a large width cig to the wavefunction of 
/j. In the limit ^ dio, dn but still ^5 «i L, A^j obtains an AFS structure 



72 \ (yg)^ 



-^\[y = -JT^f^Va^ y\) e"'^^^^ , (69) 

where is the center coordinate of the ta wavefunction relative that of the Higgs, and yl) 
is the broad- width wavefunction of centered at y-. Scanning y- and y^, the first factor becomes 
a random coefficient of order unity for all elements of the 3x3 Yukawa matrix, while the second 
factor determines the flavor suppression. As the three eigenvalues are roughly the same as 
the three flavor suppression factors, the distributions of the three eigenvalues, A^ 2,35 should be 
roughly the same as those of the three AFS suppression factors. 

Therefore, in the Gaussian landscape on the shape of the distribution functions of the 
charged lepton Yukawa eigenvalues should be like those in the left panel of Figure [61 We will 
discuss the overall range of the hierarchy shortly. If we use a common width parameter for 
10 = (g, u, e) and a different common width for 5 = (J, /)0 then the distributions of eigenvalues 
for the down quark sector should be the same as those of the charged lepton sector. On the 
other hand, the up sector Yukawa couplings involve two fields in the 10 representation, so the 
distribution of eigenvalues should be closer to those in the central panel of Figure [6] (see section [3] 
regarding the limitations of applying the AFS approximation to small eigenvalues). The results of 
a numerical simulation, displayed in Figure [161 confirm that the AFS approximation captures the 
width-parameter dependence of the distribution functions. Although we do not present numerical 
results for Gaussian landscapes other than on S^, the approximate analytic distribution functions 
for the AFS suppression factors on T^, 5*^, and are presented in Figure [TH In principle one 
can generalize the analytic discussion to many geometries. 

In the limit of broad 5 wavefunctions, ^ dio, dn (but ^ L so that values of (p\ remain 
random), the logarithmic range covered by the charged lepton Yukawa eigenvalues is 

For simplicity we choose d^o and dn to be the same, giving Alne = —{L/d)^ /16. This is 3/4 
of the result Alne = —{L/dY/12 when all of the three widths are the same. In the numerical 
results of Figure [161 the overall hierarchy in the up sector is slightly larger than that of the 
down/charged- lepton sectors, consistent with these analytic expectations. 

In Gaussian landscapes with rfg ^ c/io, dn-, \b,T tend to be predicted larger than \t. This 
is of course a problem if there is no factor such as tan (3 in the two-Higgs doublet model. This 
problem can be solved in two different ways. One is that there may be environmental selection 
in favor of a large top Yukawa coupling, which makes the observed top Yukawa coupling at the 



"'^^It is an interesting theoretical question whether higher dimensional field theories can give rise to independent 
scanning of the centers of each of the 5, u, d, I, e wavefunctions while preserving the SU(5)gut symmetric widths. 
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Figure 16: Distributions for Yukawa eigenvalues from a Gaussian landscape on S^: Qcr = 0.2 
for all overlap integrals and the Gaussian widths are d/L = 0.08 for all fields except for d and /, 
which have d^/L = 0.3. 



far upper end of its prior distribution, while the bottom and tau Yukawa couplings are not. This 
scenario does not work very well in the Gaussian landscape on S^, but it may work in Gaussian 
landscapes in higher dimensions, since the distributions of the largest eigenvalues receive a high- 
end tail in D > 1 extra dimensions (see section [5]). The other solution is that the coefficient 
Qcs of overlap integration is smaller for the down- type/charged- lepton sectors than that for the 
up-type sector. Since the overall hierarchy for the up-type sector is larger than in the down and 
charged lepton sectors, and because the distributions for the smallest eigenvalues are very broad, 
this can be done while still maintaining the order of magnitude agreement between X^^e and A^. 

It is interesting to note a difference between Gaussian landscapes with large d^ and the AFS 
models that accomplish neutrino anarchy with very weak AFS for the lepton doublets. In these 
AFS models, the logarithmic range of the hierarchy in the down- type/charged- lepton sectors is 
half of that of the up- type sector. On the other hand, this ratio is typically about 3/4 in the 
Gaussian landscape on (when d-^o = dn)- Of course the precise number 3/4 may be different 
for different numbers of dimensions, geometry, and choice of d-^o/dn- The difference between 
AFS and Gaussian landscapes stems from the fact that in the latter the AFS suppression factors 
associated with 10 are different for the up-type and down-type/charged-lepton Yukawa matrices, 

gg,n(u-sector) ^ g d{^+2d\j 2d{^ ^ gg,e(d/e-sector) ^ ^ l{d{^+d\j) _ 

In other words, the flavor suppression factors are not determined by symmetry charges intrin- 
sically assigned to fields in the 10 representation, but depend on the nature of the other fields. 
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In our universe, the down-type/charged-lepton sector hierarchy is not as small as half of that of 
the up-sector, and Gaussian landscapes have enough flexibility to accommodate this situation. 



6.3 Neutrino Mass Hierarchy 

Very small neutrino masses are obtained via the seesaw mechanism. Specifically, we assume 
left-handed neutrino masses derive from the effective dimension-five operators, 

C,f, = ^hl,hh, (72) 

which are generated after integrating out heavy right-handed neutrinos, with interactions 

C = Mc^p Vr^ Vrp + Vr^ k h ■ Cij = (A"^ A'^) . (73) 

In a Gaussian landscape the neutrino Yukawa couplings \^^^ are generated by scanning the 
center coordinates of Gaussian wavefunctions, analogous to the generation of X^^^. Yet the low 
energy observables in the neutrino sector also depend on Mcap, so we need to introduce some 
additional assumptions regarding how the Majorana mass terms of the right-handed neutrinos 
are generated. The traditional approach has been to assume an AFS pattern in the mass matrix 
Mcap. However our Gaussian landscapes are not based on the breaking of any flavor symmetry, 
but on the overlap integration of localized wavefunctions. Therefore we incorporate this principle 
into the generation of right-handed Majorana mass terms. 

Right-handed neutrinos are singlets under the Standard Model gauge group, but they are 
charged in more unified gauge groups0 For example, unified gauge groups such as SU(4)c x 
SU(2)l X SU(2)r, SO(IO), and Eg contain a B — L symmetry, and the right-handed 

neutrinos carry non-vanishing charges under this symmetry. In order for charged particles to 
have Majorana mass terms, some scalar fields on 3+1 dimensions have to be inserted to form 
gauge-invariant operators, and the expectation values of these scalar fields convert the operators 
into Majorana mass terms: 

>C = (0sb) l^Ri^R, or C= {(j)'sB(t>'sB) ^R^R- (74) 

The scalar fields (psB or (p'g^ are singlets of the Standard Model gauge groups, but carry charges 
under whatever symmetry the right-handed neutrinos are charged. In other words, they are 
moduli fields describing the symmetry breaking of the more unified symmetry. 

Since the Majorana mass terms fTMl) involve vacuum expectation values of scalar fields, the 
coefficients of the mass terms will involve overlap integrations: 

Mc„^(x / dy^'^^^l^ipy, (75) 



-'^^Tlie term "unified gauge group" is here used in a loose sense. For instance, we do not assume that the 
Standard Model gauge group is unified as an effective field theory in 3+1 dimensions. 
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where Lp^J^ is the zero-mode wavefunction of the a copy of the right-handed neutrino and tp 
represents all of the other effects of localization in the extra dimension (s), including the wave- 
functions of the symmetry-breaking moduli fields (j)sB and/or (p'g^- These fields have nothing to 
do with the Higgs boson, and hence there is no reason to believe that ^p^^ is localized at the same 
place in the extra dimensions as the Higgs boson wavefunction. Thus while the AFS structure of 
all of the Yukawa couplings A"''''^''^ is due to an overlap between fermion wavefunctions and that 
of the Higgs, any flavor structure in the Majorana mass term follows from an overlap involving 
ip^^ . In particular, unless there is a strong correlation between the wavefunctions of ip^^ and 
the Higgs, the flavor structure of the mass matrix Mcap is statistically independent of that of 
the Yukawa matrices. 

We study the consequences of this new flavor structure in sections 16.3.11 and 16.3.21 However, 
before we begin some comments are in order. Our first remark concerns the overall mass scale 
M of the right-handed neutrinos. The overlap integration in fl75|) sets the flavor structure of 
Mcapi but does not say anything about the overall scale of the symmetry breaking. Thus 
an extra assumption has to be introduced in order to set the distribution of the overall scale 
M. Even if we know this distribution from theoretical considerations, it may be modulated by 
environmental selection related to leptogenesis. Therefore we set aside this issue and instead 
focus only on the hierarchy among neutrino masses and mixing angles in the lepton sector, which 
can be determined independent of the overall scale of the neutrino masses. This is equivalent 
to studying a fixed M cross section of the full landscape. Our second comment concerns the 
form of the wavefunctions of the right-handed neutrinos and the symmetry breaking source ip^^ . 
To date we have very little knowledge about these fields; yet in a Gaussian landscape their 
most crucial aspect will be whether in the extra dimensions their wavefunctions are localized or 
not. Therefore we set aside the theoretical origin of these fields and simply represent them with 
Gaussian wavefunctions fl2Ul) . while considering the width parameters and as unknown. 
By choosing dj^i and to be large or small, we can simulate various possibilities. For example, 
Majorana mass terms generated by world-sheet instantons wrapped on topological cycles may 
be mimicked by a (possibly multi-centered) Gaussian wavefunction with a small dsB- 

It turns out that the mass matrices of the neutrino sector are quite different depending on 
whether the right-handed neutrinos have localized wavefunctions or not. Therefore sections 16.3. II 
and l6.3.2l are separately devoted to these two possibilities. In either case, to generate a statistical 
ensemble of Ajj^j and Cap the center coordinates of wavefunctions are scanned randomly and 
independently for ip^^ , the right-handed neutrinos, lepton doublets, and the Higgs. 
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6.3.1 Narrow Right-Handed Neutrino Wavefunctions 

If right-handed neutrinos have narrow wavefunctions, then the neutrino Yukawa matrix has a 
structure similar to that of charged leptons. That is, we have 

This is f l69|) with d\o replaced by c/at; and again the Qai are effectively random coefficients of 
order unity. The range of the hierarchy is Alne'^ = — L^/[8((i^ + c?^)]. 

Meanwhile, the Majorana mass matrix c„/3 of right-handed neutrinos is almost diagonal, 
because the off-diagonal entries are suppressed due to the (statistically) small overlap of wave- 
functions between different right-handed neutrinos. If the symmetry-breaking source ip^^ is also 
very localized, then the Majorana mass matrix has the approximate structure 

Cap ~ 5ap e ''sB+-^N = , (77) 

where cise is the width of Gaussian wavefunction of ^p^^ and y^^ is the center coordinate of this 
wavefunction. The positions y^^ and y^^ are scanned randomly and the largest, middle and 
smallest become the three eigenvalues of the right-handed neutrino mass matrix. In particular, 
their distribution functions should be like those of the AFS suppression factors ei 2,3 in Figure [6] 
or [T5l For Gaussian landscapes on we have Alne^^ = — L^/[4(2(i|^ -|- d\i)]. Results of a 
numerical simulation on the landscape are shown in the first row of Figure [T71 confirming the 
theoretical expectations so far. If the symmetry-breaking source ip^^ does not have a localized 
wavefunction, then we can take dsB to be very large. Then Alne*^ ^ 0, meaning that there is 
not much hierarchy among the mass eigenvalues of right-handed neutrinos. 

From fl75]) . fITBl) . and fl77|l the mass matrix of low-energy neutrinos is approximately 

Cij ~ X] gaigaj-J^ ■ (78) 
a " 

Since we introduce no correlation between the Higgs and (p^^, we do not expect net statistical 
cancellations in the ratio {^aY/^a- This contrasts sharply with flavor theories that assume an 
AFS structure for the Majorana mass matrix of right-handed neutrinos. Let the distributions of 
In and In be denoted by 

dPiy"") = foiy'') dy^ , dP{y'') = fM{y'') dy^ , (79) 

where we introduce the new variables y-^'^ = Ine^'^'^/Alne*^ {y^ is normalized to run from zero 
to one; note that the variable y*^ does not necessarily run from zero to one). Because of the form 
of and e*'^ in (1751) and fl77|) . foiv) and fuiv) are given by the volume distribution function 
f{y) of (l6Up . after re-normalizing the distribution function and rescaling the argument, if nec- 
essary. To approximate the distributions of eigenvalues of In Cij, we first introduce the variable 
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logio(mi/m^) 



logio(m2/m^ 



Figure 17: Distributions of the eigenvalues of the right-handed neutrino mass matrix and the 
left-handed seesaw Majorana mass matrix, normalized by M and m^, = (h)'^ /M, respectively. 
The Gaussian landscape on is used for this simulation, with d/L = 0.08 for the wavefunctions 
of T^Ra, 1 and the Higgs while d/L = 0.3 for the lepton doublets. 



z^^ = In [(e'^)^/e^] /Alne^. The distribution of follows from appropriately convoluting the 
distributions of Ine^ and Ine*^, 

^^P- = /ss(^^^) = llfD (^/2) fuiz - n dz . (80) 

Note that the variable z^^ runs from —2{d'jj + (i^)/(2(i|^ + dj^) to 2 on S^. The distribution 
function fgg can be calculated explicitly once the underlying volume distribution function f{y) 
is known; the results for the Gaussian landscapes on and are shown in Figure [HI Using 
the volume distribution function fssi^ss), we obtain the distributions of the three eigenvalues of 
Cij by integrating the other two variables out of the combined distributior0 

dP{zl^ zf, zf) = 3\UzT)Uzf)fss{zf) e{zf - zf) Q{zf - ^r) dzl^dzl^dzf , (81) 

where ^^^23 ~ ~8[((i^ + (i^)/L^] ln(mi^2,3/"^!y); with m,^ = (h)'^ /M. Distribution functions of 
other observables such as a ratio of low-energy neutrino masses are also calculated from this 
combined probability distribution. Some examples are displayed in Figure [181 Note that the 
low-energy neutrino masses tend to have a very large hierarchy, because the distribution function 
fss covers the enormous logarithmic range 

A,„™„ ^ 2A„„- . Al,„« ^ -1 (-^ . . (82) 



"'^^The combined probability distribution of neutrino mass eigenvalues is discussed in |10| for the case with 
fo and /a/ proportional to a (5-function. Reference [13j deals with the neutrino anarchy scenario, where the 
diagonalization effects on the eigenvalue distribution are also taken into account. 
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Figure 18: Distributions of density functions fss{z^^) (left column), low-energy left handed neu- 
trino eigenvalues (middle column), and the mass ratios (mi/m2) and {m2/m^) (right column). 
The top row is for while the bottom row is for S*^ (to reduce computation time we have 
approximated with a polynomial that is always within a part in a thousand of the corre- 
sponding expression in (!6T|) ). In the top row the distributions for the middle eigenvalue m^jvny 
and the ratio m2/m^ are displayed with a dashed line. The two mass ratios are the same on S*^. 
For these curves du/L = dpf/L = dss/L = 0.08. 



Specifically, the hierarchies of neutrino Yukawa matrix and the Majorana mass matrix add 
rather than cancel, c.f. [HI [13] • Even when dsB is very large, low-energy neutrino masses have a 
hierarchy twice as large as that of eigenvalues of the neutrino Yukawa matrix. 

The low-energy neutrino eigenvalues of Figure [17] show a remarkable similaritjQ in shape with 
the corresponding distributions in Figure [TBI confirming that the theoretical derivation of dis- 
tribution functions is fairly reliable. Although results of a numerical simulation on the Gaussian 
landscape on are not presented here, we have confirmed that the distributions of \n{mi/my) 
are similar to corresponding distributions in Figure [TSl This success arises because the overlap 
integration for A^j, AJ^^ and c^/? involves only two very localized wavefunctions. Although we 
have seen around fl37p that the analytic distribution functions based on the volume distribution 
function have a systematic error due to the compactness of the extra dimension(s), the error 
occurs only for overlap integrations involving three localized Gaussian wavefunctions. Therefore 
these analytic results can be used to compute distribution functions of flavor observables in the 
lepton sector in other geometries for which numerical simulation is costly. 

The underlying geometry is reflected in the volume distribution functions fn and /a/. How- 

-'^^The peak of the distribution of logj^Q m2/m^ in Figure [17] is near —2, as opposed to the expected location 0, 
shown in Figure [TSl The peak of m2/m^ being lower than that of ms/rrii, is not surprising. We have seen similar 
phenomena in Figures [5l [TTl [Til and [TBI As mentioned earlier, these are likely due to the diagonalization effect. 
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ever, these are integrated once to obtain /ss, which smears the effects of geometry. The original 
/'s are integrated at least once in the distribution functions of ln(mj/m,y), like those of the 
up-sector eigenvalues. Moreover, the distribution functions of observables such as ln(mi/m2) 
and In(m2/m3) involve two integrations over the /'s, such that very little geometry dependence 
is left in the distribution of the observable!^ 

6.3.2 Broad Right-Hand Neutrino Wavefunctions 

If the zero-mode wavefunctions of right-handed neutrinos are not localized, i.e. the width pa- 
rameter (In is not much less than L, then in the neutrino Yukawa coupling only the Higgs has 
a localized wavefunction. Thus the overlap integral is evaluated around the peak of the Higgs 
wavefunction. Expanding the other wavefunctions. 



where all of the wavefunctions are evaluated at y . Note that the first term is a rank one matrix 
while the additional terms are suppressed by (dj^/minjdAr, cig})^. Thus we obtain a neutrino 
Yukawa matrix that has the form of the "democratic" mass matrix ansatz, i.e. it is rank one at 
leading order [18]. Deviations from absolute democracy (i.e. from being rank one) result from 
the wavefunctions for U and Vj^a not being absolutely flat. The first row of Figure [19] shows the 
distributions of the three eigenvalues of the neutrino Yukawa matrix. Note that y?, ^p' and if" 
have upper limits and there are no effects enhancing the middle eigenvalue; hence the sharp cut- 
off of the distribution of In Ag. An interesting conclusion is that the neutrino Yukawa matrix has 
a certain amount of hierarchy, no matter how large or small the width parameter oi the right- 
handed neutrinos is. Small Yukawa couplings result from small overlap integrals for localized 
right-handed neutrinos, or from rank reduction for non-localized right-handed neutrinos. 

The Majorana mass matrix of right-handed neutrinos is not hierarchical when the source of 
symmetry breaking ip^^ is not localized. Thus in the second row of Figure [19] the distribution of 
M2/M is roughly order unity, and the hierarchy M2/M3 is not more than what is expected from 
the d,ago„al,;ation effect§ On the other hand, the Majorana masses of right-handed neutrinos 

-'^^The cuspy peaks in the distributions of fss and ^^{mi/m^) on are only logarithmic; the y^^^"^ singularity 
of the original /51 is integrated at least once in each of these. We chose a fine binning in Figure [17] in order 
to accentuate the logarithmic singularity and emphasize the agreement between the simulation and analytic 
derivation. When a coarser binning is chosen, the logarithmic singularities are smeared out. Thus, there is not 
much practical difference in the distributions of ln(mi/m,y)'s whether or or any other geometry is used. 

^"We do not have a clear explanation why the distribution of Mi/M has a long tail toward smaller values. 






we see that the neutrino Yukawa matrix is of the form 
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logio(A/i/M) logio(M2/M) logio(M3/M) 
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Figure 19: Distributions of the eigenvalues of the left-handed neutrino Yukawa matrix, the 
right-handed Majorana mass matrix, and the low-energy left-handed Majorana mass matrix. 
The numerical simulation uses a Gaussian landscape on S^, with dfj/L = 0.08 and dss/L = 
dj^/L = d^/L = 0.3. 

become hierarchical when ip^^ is localized. Then Cajj is rank one at leading order, and deviation 
from being rank one should come at the order of {dss/ d^Y . 

The third row of Figure [12] shows the distribution of mass eigenvalues of the low-energy seesaw 
neutrinos for d^jL = dss/ L = 0.3. The distribution of ln(m2/m,^) is sharply cut off at its largest 
value; this cut-off is presumably traced back to the similarly sharp cut-off of the distribution 
of the second largest eigenvalue of the neutrino Yukawa matrix Ag. As opposed to when the 
right-handed neutrinos have localized wavef unctions, it is not as straightforward to develop 
an analytical understanding of the distribution of seesaw neutrino mass eigenvalues when the 
right-handed neutrinos have non-localized wavefunctions. Instead, we run numerical simulations 
with the Gaussian landscape on S^, using a range of dj^/L G [0.2,0.5] and dss/L G [0.08,0.5]. 
We find that the distributions of ln(mj/mi,) remain qualitatively the same as in the last row of 
Figure [TUl It is only when dsn is much smaller than d^ that there is a qualitative difference; then 
the Majorana mass matrix of right-handed neutrinos has hierarchical eigenvalues, and lighter 
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right-handed neutrinos give rise to enhanced values of m2/m^, erasing the sharp cut-off in the 
distribution of \n{'m2/m^). 

6.3.3 Neutrino Mass Hierarchy in Gaussian Landscapes 

It is a robust consequence of Gaussian landscapes on any geometry that the seesaw masses 
are hierarchical. If the right-handed neutrinos have localized wavefunctions, then the neutrino 
Yukawa matrix obtains statistically hierarchical mass eigenvalues just like in the analysis of 
sections [2] and [31 On the other hand, if the wavefunctions are not particularly localized over the 
extra dimension(s), then the Yukawa matrix always has small eigenvalues due to rank reduction. 
Because of the absence of correlation between the Majorana mass terms of right-handed neutrinos 
and the neutrino Yukawa couplings, the hierarchical structure of the neutrino Yukawa matrix 
remains in the seesaw mass eigenvalues. As a consequence, an inverted hierarchy is unlikely. In 
addition, it is likely that nii <^ m2 and 



Although Gaussian landscapes predict that the seesaw masses are hierarchical on average, 
the distribution of mass ratios from Gaussian landscapes are qualitatively different, depending 
on the choice of the width parameters and dsn- Consider the distribution of m2/m^ for 
various choices of these width parameters, displayed in Figure [201 When neither the right- 
handed neutrino nor the symmetry-breaking wavefunction(s) are localized, the distribution of 
ra^jm-i is sharply cut off at around 10~^, seen in Figure [201 (d) and (e). In this case there is no 
chance to accommodate the observed value. 



If (ijv is large and the wavefunction i^^^ is strongly localized in the extra dimension(s), then 
the distribution of m2/m^ has a tail extending upward as in Figure [201 (f). This is because of the 
previously mentioned enhancement in the seesaw mass matrix — due to small eigenvalues in the 
Majorana mass terms of the right-handed neutrinos — which erases the sharp upper bound on X^. 
The distribution in Figure [201(f) barely covers the observed value ( |87l) . but the width parameters 
may be chosen a little differently to provide better coverage in the tail of the distribution. Yet 
even in this scenario, the distribution of m2/m^ is peaked at a value much smaller than the ob- 
served value ([871) . If environmental selection factors (such as those associated with leptogenesis) 
and weight factors associated with cosmological evolution are in favor of such an outlier, then 
this scenario may be acceptable. Otherwise, one should conclude that this scenario is not a good 
approximation to the actual landscape governing the flavor structure of the multiverse. 

It is left to consider localized wavefunctions for the right-handed neutrinos. The approximate 
distribution of m2/m^ is derived analytically in this scenario, and plotted in the right column of 
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Figure 20: Distributions of log^^gl'^s/'^s) for various choices of parameters. The Gaussian 
landscape on 5*^ was used for all simulations except (g), for which was used. In (f), dss is 
chosen to be very small so that djj/d'^g^ is large, while (h) involves complex Gaussian profiles 
(see section I6l4|) . Unless otherwise specified, d\o/L = dn/L = 0.08 and d^/L = 0.3. 



Figure [181 Regardless of the geometry this distribution covers a broad logarithmic range, with 
the greatest weight toward the largest value, 7722/7723 ~ 1. These analytical distributions capture 
the qualitative features of Figure [20] (a), (b), and (g) very well, with one important difference: 
the distributions in the numerical results fall off to almost zero for m2/m^ ^ 10~^-10~^. This 
behavior of the numerical results is presumably due to the diagonalization effect, a phenomenon 
first mentioned in section 13.21 and subsequently referred to with regard to Figures [3, [TT], [TU [TBI 
and[T71 This reduced weight in the probability distribution functions at 7712/7773 ^ 10"^ means 
that this scenario is also in conflict with the observation flHTl) . 

There is an important caveat, however, which is that the diagonalization effect changes 
as the real-valued matrices that we have dealt with so far are generalized to complex-valued 
matrices. See the appendix of [13], for example. Figure [2D] (h) is based on an extended version 
of the Gaussian landscape on (presented in section 16.41) that includes complex phases. The 
diagonalization effect clearly has much less impact on the distribution function of 7772/7773, and 
this distribution extends almost all the way up to 7712/7773 ^ 1. Thus we find the scenario 
with localized wavefunctions for the right-handed neutrinos is compatible with the large value 
observed for IS.mf^/ ^ml^^, assuming the Gaussian landscape contains complex phases. 
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Figure 21: The distributions of the PMNS mixing angle 6*23, the (2,3) element of the diagonal- 
ization matrix V^, and that element of V^. The distributions result from a numerical simulation 
of the Gaussian landscape on 5*^ with d/L = 0.08 for all fields except d^/L = 0.3. 




^23 ^23 

Figure 22: Distribution of the (2, 3) component of the diagonalization matrices and in 
the complex plane. The simulation uses the Gaussian landscape on 5^, with d/L = 0.08 for all 
fields except rfg/L = 0.3; while r = 3. The distributions of jV^gl and IV23I remain much the same 
as in Figure [21], but the phase angle prevents precise cancellations in ^23- 



6.4 Mixing Angles and CP Phases 

We introduced non-localized wavefunctions for the fields in 5 because we expected this to result 
in large leptonic mixing angles. However numerical simulation reveals that the mixing angles, 
particularly ^23? are still very small, as can be seen in the left panel of Figure HTJ The figure 
presented here is a simulation using small d]^/L and dsB/L; but the 6^23 distribution is confined 
near 6'23 ~ no matter how small or large these parameters are taken to be. 

The distributions of V23 and V23 in Figure [21] explain why, contrary to our expectation, we 
failed to get large 6^23. An SO (3) matrix acts on the left-handed charged leptons e^j C li 
when diagonalizing the charged lepton Yukawa matrix A^j, and acts on the left-handed 
neutrinos un C k when diagonalizing the seesaw mass matrix Cij/M. Their (2,3) components 
are large, just as we expected. Yet the mixing angle of atmospheric neutrino oscillation, the (2,3) 
component of f/ = {V'^)~^V'^, is small due to a precise cancellation between the large V23 and 
V23. This cancellation is not a new phenomenon. The democratic mass matrix ansatz, originally 
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applied to the quark sector, yields small mixing angles through this precise cancellation. 

To obtain large mixing angles, some modification is needed to prevent this cancellation. We 
replace the real-positive valued Gaussian wavefunctions in (120]) by Gaussian wavefunctions with 
a winding complex phase @ 



introducing a new parameter r. For simplicity we use the same value of r for all of the wave- 



be the complex conjugate of that of h, since for example in the framework of section [7] these 
wavefunctions are not just scalars but have more complicated internal structure!^ 

As seen in Figure [22], using this form of wavefunction makes and V23 complex, and their 
phases decrease the likelihood of a precise cancellation in 623- Furthermore, these complex Gaus- 
sian wavefunctions introduce CP violation into the fiavor physics. Meanwhile, this introduction 
of complex phases has little effect on the distributions of charged fermion mass eigenvalues. 
This can be understood analytically by tracing the discussion in the preceding sections with 
the wavefunction (l88l) . We also confirmed this result by running numerical simulations. The 
diagonalization effect is not taken into account in the theoretical arguments, but in practice this 
matters only for the m2/m^ neutrino mass distribution, and we have seen that this improves the 
agreement with observation. 

On the other hand, introducing complex phases dramatically changes the distributions of 
mixing angles. Figure [23] shows the distributions of the three leptonic mixing angles resulting 
from a numerical simulation. We see that the mixing angles of solar and atmospheric neutrino 
oscillations, 612 and ^23, are likely to be any value between and 7r/2. Distributions of the 
oscillation amplitudes sm{29ij) are almost flat for both of these mixing angles. Peaks at maximal 
mixing come from the Jacobian [13]. Comparing the distributions of ^23 between Figure [23] and 
Figure [2T] we see that the introduction of complex phases (and hence CP violation) is an essential 
ingredient in obtaining large (and even maximal) mixing angles in neutrino oscillations. 

The mixing angle ^13 can also be of order unity, which is similar to the situation in neutrino 
anarchy [9] [13]. In the U(3)-invariant measure to be expected in neutrino anarchy, however, the 
distribution function of sin(26'i3) is peaked at sin(26'i3) = 1, and vanishes at sin(2^i3) = 0. Thus 
all three mixing angles in the lepton sector are likely to be large. In particular, according to 
neutrino anarchy using the U(3)-invariant measure, about 60-90% of the statistics are already 

^^This form for the wavefunction is not without motivation, c.f. the discussion surrounding (|154p . However, we 
use this wavefunction only as one of the simplest means to introduce complex phases into Gaussian landscapes. 
There are more complicated, more sophisticated, and possibly better justified ways to introduce complex phases 
into Gaussian landscapes, or landscapes based on extra dimensions in general. It would be interesting to explore 
such models that include complex phases to study the correlation between the induced CP violation and the 
distribution of other obscrvablcs. However this subject is beyond the scope of this paper. We restrict attention 
to (|88[) and find that even this simple model is sufficient to obtain large leptonic mixing angles, when df, is large. 

^^We have checked that one can obtain results very similar to those presented in this section, but taking 
ip^ — i^ip^ )*, if one also allows for non- universal values of r. 





functions in the overlap integrations, tp' 



,q,u,d,i,e,u,SB,h,h* _ rjj^g wavefuuctiou of h* does not have to 



51 




Iogio(2^i2/7r) 



logio (26^23 vr) 



logio sin 6^13 



Figure 23: Distribution of the leptonic mixing angles in the Gaussian landscape on S^, with 
complex-valued wavefunctions with r = 3 and d/L = 0.08 for all fields except d^/L = 0.3. 



excluded by the current experimental limit on ^13 [HI [13] • In the Gaussian landscape on 5"^, 
on the other hand, the ^13 distribution is weighted more toward sin(26'i3) = 0. As we see in 
section 16. 5[ only about 30% of the statistics have been excluded. 

For comparison, the distributions of the quark mixing angles are also simulated in the same 
Gaussian landscape and are presented in Figure [241 We see that the complex- valued Gaussian 
profiles maintain small mixing angles when the overlap integrals involve narrow Gaussian widths 
for lO's and the Higgs. Interestingly, simply choosing different widths for particles in the 5 and 
10 representations, in the presence of large CP-violating phases, allows for very different flavor 
structures between the quark and lepton sectors. Thus pursuing a microscopic description behind 
the statistical distributions of Yukawa couplings has enabled us to go beyond the results of p[TO]. 

The distributions of the CP violating phases are also shown in Figure [251 Our parametrization 
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Figure 24: Distribution of the quark mixing angles in the Gaussian landscape on S*^, with 
complex-valued wavefunctions with r = 3.0, dio/L = du/L = 0.08 and rfg/L = 0.3. 



of the CKM and leptonic mixing matrices V and U is quite standard: 

C12C13 S12C13 «,.p-i<5cKM 

V = I -S12C23 - Ci2S23Si3e^^CKM C12C23 - Si2S23Si3e^^CKM ^^^^3 | , (89) 

S12S23 - Ci2C23Si3e*'^CKM -C12S23 - Si2C23Si3e^''cKM 




C12C13 S12C13 
U = I -S12C23 - Ci2S23Si3e^'''' C12C23 - si2S23Si3e^'''^ S23C13 | | e^"2/2 I ^ (90) 

■S12S23 ~ Ci2C23Si3e*'^'' — C12S23 — S 120235136*'''' 

where Sij = sin Oij and = cos 6ij. The Majorana mass phases of the neutrinos, ai,2 = 
Arg(mi, 2/^^13), have almost flat distributions. On the other hand, the distributions of the two 
other CP phases, 5ckm and 6^, have peaks at 5 ~ on top of otherwise flat distributions. We 
have not studied where this structure comes from. The scatter plot of V^g in Figure [22] reveals 
some structure — the scatter is dense in some places and thin in others — and the peaks at 5 = 
in the CP phase distributions may have something to do with this structure, which presumably 
originates from the specific form of the wavefunction (188|) . If this guess is correct, then the peak 
plus flat structure in the CP phase distributions is likely to be an artifact of the particular way 
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Figure 25: Distribution of the CP phases in the quark and lepton sectors, based on the Gaussian 
landscape on 5*^ with r = 3.0, d/L = 0.08 for all fields except d^/L = 0.3. 



in which we have introduced complex phases into this landscape. It would be interesting to see 
how the CP phase distribution changes when complex phases are introduced in a different ways, 
but this subject is beyond the scope of this paper. Although we expect that the peak structure 
may deform, disappear, or become less significant, we consider it unlikely that the flat part of 
distribution would disappear and that CP phases would be predicted to be very small. If this 
expectation is correct, then these landscapes are consistent with the measured value of the CP 
phase in the quark sector. 

6.5 Conditional Probabilities 

The neutrino sector has three observables yet to be measured, sin ^13, 61, and rupp. It would 
be interesting if landscape approaches to understanding flavor could make predictions for these 
observables. We have developed our Gaussian landscapes so that all the observables measured 
in this universe are not too atypical. In this section, we use the Gaussian landscape on S^, with 
the set of parameters r = 3, dfj/L = d^o/L = djy/L = dss/L = 0.08, and d^/L = 0.3, as an 
initial example to explore what kind of predictions can be obtained from this approacho 

The mixing angle ^13 is typically of order unity, and there is no reason to expect from this 

^^This Gaussian landscape has 19 scanning parameters, corresponding to the center coordinates of the wave- 
functions of the quarks, leptons, Higgs boson, and symmetry-breaking field(s) tp^^. Meanwhile, this landscape 
now predicts distributions for all 22 flavor parameters. Although this implies that the landscape makes three pre- 
cise (zero-width) predictions; this is a result of our arbitrary choice to consider parameters such as g, i/dg.io.Af.SBi 
and r as fixed. Specifically, these predictions disappear when these parameters also have probability distributions 
around the values we fixed. Therefore we do not pursue the possibility of a truly predictive Gaussian landscape. 
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Gaussian landscape that it is very small. Indeed, the current experimental limit has excluded a 
significant fraction of the ensemble of vacua: 



Pfsin^i3 > 0.18) 



33%. 



(91) 



Nevertheless, the experimental limit is not very strong and so a significant fraction of vacua — 
67% of the total ensemble in this Gaussian landscape — sits within the experimental bound. On 
the other hand, future experiments are expected to have a sensitivity down to about 6^13 ~ 10"^. 
The probability that ^13 is too small to be measured by such future experiments is rather low 
within this Gaussian landscape; we find 



P(sinei3 < 10 



12%. 



(92) 



The second yet-to-be measured parameter is the CP phase of neutrino oscillations. The 
probability distribution of this is displayed in the bottom left panel of Figure [251 Finally, 
neutrinoless double beta decay measures |«^^^|, where 



AR = ^ 



i=l 
Uel 

The second term of (l93l) is roughly 

Ue3 



■.rui 



171211^2 



Ue2j fn2 



nil 

7712 



(93) 
(94) 



e2 



7713 



m2 



0.57 X 



sin 6' 



13 



0.18 



0.18 X 



sin 6 



13 



0.10 



(95) 



after using (l86ll and (l87j) . Once sin 6^13 is measured, the second term of (p3|) can be estimated. 
Yet we still need to know the last term, AR, to make a prediction for |?7i/3/3|. The Gaussian 
landscape distribution of \AR\ is found in the top right panel of Figure [261 

We can ask more specific questions to further constrain the probability distributions of these 
three observables in universes like ours. That is, when making predictions for future measure- 
ments, we can condition the probability distributions based on quantities that we have already 
measured, as opposed to using the a priori probability distributions in the landscape. However, 
if we were to impose the current experimental limit on all lepton fiavor parameters, an enormous 
numerical simulation would be required, since most of the simulated universes would not pass the 
cut. This does not imply that our universe is very atypical, only that the present experimental 
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BnCnD 




61% 


75% 


78% 


77% 


94-95% 


79% 


87% 


96 ± 1% 



Table 2: Probability to measure P(sin^i3 > 3 x 10~^) under various cuts. Up to ±1% uncertain- 
ties arise from the limited statistics gathered in the numerical simulation. We use the Gaussian 
landscape on with r = 3 and d/L = 0.08 for all fields except d^/L = 0.3. 



error bars are small. Therefore, for practical reasons we employ much looser experimental cutJ^ 



A 
B 
C 
D 



10^' < Ami/Ami^ < 10^\ (96) 

sin2(2^i2) > 0.7, (97) 

sin2(2^23) > 0.8, (98) 

sin0i3<O.18, (99) 



and study whether these cuts influence the distributions of 6*13, AR and Si^. 

Table[2]shows how the probability to measure sin ^13 > 3x10^^ changes depending on whether 
or not the cut conditions A-C are imposed. The loose cuts A-C have a significant impact on 
the prediction. Among all the samples that are consistent with the current experimental limit 
on sin ^13, 40% of them are below 3 x 10"^. However, after imposing various other experimental 
constraints, A-C, only a few percent have sin6'i3 < 3 x 10~^. Therefore, for the wavefunction 
parameter choices used in this section, the Gaussian landscape on predicts a very high 
probability of measuring 6*13 in future experiments in our universe. 

When sin ^13 is large, there is a good chance to observe CP violation in neutrino oscillation 
experiments. The a priori vacuum statistics for this parameter are peaked at 6^ = (bottom 
left panel of Figure [25]) . However the distribution of 6^, changes as some of the cuts A-D are 
imposed. The impact of these changes is clearly visible when the mixing angles of both the solar 
and atmospheric neutrino oscillations are required to be large; then the peak of the CP phase 
distribution is significantly reduced (bottom middle panel of Figure 1251) . This is reasonable, 
since complex phases play a crucial role in avoiding cancellation between large mixing angles in 

and V^. In scatter plots in the sin^(2^i3)-(5^ plane (Figure [26!) . we see that the distribution 
of the CP phase becomes uniform after cuts A-D are imposed. Therefore, after conditioning 
distributions based on our loose "experimental cuts," both sin ^13 and the CP phase become 
more likely to be larger, and hence more likely to be discovered in future experiments. 

The loose cuts A-D also have some impact on the distribution of |Ai?| (Figure The 
condition B alone, requiring \Ue2\ be large, removes almost all the distribution log^^Q |Ai?| > 



^^One might consider these cuts as a tool to study the effects of possible cosmological or environmental selection 
in the multiverse, just like a large top- Yukawa cut (|43p was employed in section [3] in an attempt to study the 
impact of environmental selection in favor of a large top Yukawa coupling. Here, however, we know of no evidence 
that the weight factors are in favor of the range of parameters selected by the cut conditions A-D. 
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Figure 26: Distributions of observables in the neutrino sector for different subsets of the total 
ensemble. The Gaussian landscape on 5"^ is used for this simulation, with r = 3 and d/L = 0.08 
for all fields except d^/L = 0.3. In the log^o sin^(26'i3)-5i, scatter plots, the same number of 
points are displayed in all three figures. 



as well. Examining the samples that pass all cuts A-D, we find that 



P(\AR\ > 1] 
P(\AR\ > OA] 



(0.3 ±0.2)%, 
(9 ± 1)%, 



(100) 
(101) 



where the uncertainty comes from limited statistics in the numerical simulation. Therefore, it is 
very unlikely that AR is so large as to be significant next to the first term m2f/g2 in mp^. Hence 
this landscape predicts m^^ ^ '^2f^e2 + "^s^es ^nd 



e2l 



0.18 



sin 6 



13 



0.10 



e2l 



X 



1 + 0.18 



sin 6 



13 



0.10 



(102) 



where \m2U^ 



e2l 



(2-3) X 10-3 eV. 

We have seen that the loose cuts A-D have a significant impact on the distributions of 
sin6'i3, 6i, and |Ai?|. This study gives us a feeling for how much landscape predictions can 
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change when various weight factors multiply the simple vacuum statistics, or when various 
experimental measurements are used to condition predictions based on what we know about 
our universe. Details about these predictions will depend on the specific weight factors or how 
severely the experimental cuts are imposed. However, as we have seen, the impact of these cuts 
can be understood qualitatively, and we expect similar results for similar weights or cuts. 

Although the width parameters of 10, 5, and the Higgs have been optimized to fit the 
observed data, there is still much room left to change and dsB (or even r or how complex 
phases are introduced into the Gaussian landscape). If dN and/or dsB were chosen slightly 
larger, then the distribution of 7712/7723 covers a smaller range, and the cut A would have less 
impact while the other three cuts on the mixing angles would remain as important. Since it 
was cuts B n C and B, respectively, that had impacts on the distributions of 6^, and |Ai?|, the 
"predictions" on 61, and \AR\ would not be affected very much. On the other hand, since the 
AnB cut was important in dragging the distribution of sin 6*13 upward, the distribution of sin6'i3 



may shift downward for a slightly larger d^ or dsB- T 
be interpreted carefully, with these and other caveats 



lus the "predictions" in this section should 
in mind. 



7 Yukawa Couplings from Super Yang— Mills Interactions 

In the preceding sections, we introduced a number of toy models, collectively termed Gaussian 
landscapes, to generate statistical ensembles of Yukawa matrices. The pairing and generation 
structures of the CKM matrix along with the large mixing angles of the lepton sector were 
successfully explained as random selections from these ensembles. The most crucial ingredient 
to these Gaussian landscapes was that Yukawa couplings are given by overlap integration of 
localized zero-mode wavefunctions on extra dimensions. The correlation necessary to account 
for generation structure arose because the wavefunctions of quark doublets and the Higgs are 
relevant for both the up-type and down-type Yukawa matrices. We also introduced a number of 
additional assumptions concerning which parameters are scanned and which not. In this section 
we give give greater theoretical consideration to these assumptions. 

Consider a supersymmetric Yang-Mills theory on a higher dimensional manifold. The inter- 
nal space of D extra dimensions is denoted X^. The gauge group G is chosen so that it contains 
a subgroup H such as SU(3)c x SU(2)l x U(l)y of the Standard Model or a unified gauge group 
such as SU(5)gut- The gauge fields Am of the theory may have non-trivial configurations on 

^^For example, landscapes of Yukawa couplings originating from super Yang-Mills theory in higher dimensions 
have basis-independent distributions, but Gaussian landscapes — simplified versions of the former — do not (see 
section [7]) . This situation is compared to the basis independent implementation of neutrino anarchy in |13j and 
the original implementation in [5]. The latter predicted anti-correlations between mixing angles in the lepton 
sector, but these angles are uncorrelated in the former. Thus some of the correlations in Gaussian landscapes 
may be spurious. Nevertheless, the disappearance of correlations in the basis-independent neutrino anarchy was 
related to invariance under the U(3) group. Unless there is a similar underlying symmetry, we expect that the 
probability distributions of Gaussian landscapes are not qualitatively different from those of higher dimensional 
gauge theories. 
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Xd as long as they are stable (and hence satisfy the equations of motion). Such a gauge field 
configuration can break the underlying gauge symmetry G down to H. 

In this context the fields {fj, = 0,1,2,3) in f) become the gauge fields of the Standard 
Model or the unified theory. Meanwhile, the Kaluza-Klein spectrum of Am [m = 5, ■ ■ ■ ,4 + D) 
in g, under the gauge field background, may have zero modes (i.e., massless modes), and such 
quantum fluctuations become scalar fields in the four- dimensional low-energy effective theory. 
Gauginos in g may also have zero modes; since a spinor in a higher-dimensional space is bi- 
spinor of the Lorentz group SO (3, 1) and SO(D), and the zero-modes of the Dirac equation on 
Xd become massless fermions in the effective theory. These low-energy scalar and spinor fields 
may be charged under the Standard Model or unified theory gauge group, or they may be singlets 
under these gauge symmetries. This is worked out by decomposing the adjoint representation 
of Q into pieces irreducible under H. Charged fields arising in this way may well be identified 
with quarks, leptons and the Higgs. 

Fermion and scalar fields arising from the super Yang-Mills fields in higher dimensions have 
Yukawa couplings in the low-energy effective theory. To see this more explicitly, consider the 
Lagrangian of a super Yang-Mills theory in higher dimensions: 



r - 1 



MP _ ^„.„ MP _ M,^ 



tr {F^^F^ + 2^ tr (F^^F'^'") + ^ tr (*r^9^^) + 



9* ' 9: ' 9i 



(103) 



These three terms become the kinetic terms for the vector, scalar, and spinor fields in the effective 
theory. Meanwhile, the gauge coupling constant of the low-energy effective theory is 

1 mPv 



9eS. 9l 



(104) 



where V is the volume of the internal space. When the zero- mode wave functions <^m{y) and 
^^{y) are normalized such that 



MP d''y\ipm? = l, MP d^y|^/;|2 = l, (105) 

then canonically normalized kinetic terms result when the low-energy fields are related to those 
in the higher dimensional theory by 

Ayn{,x,y) = g^ipm{y)<p{x) , {x, y) = g^ij{y) ■ x{x) . (106) 

Then low-energy Yukawa couplings originate from the gauge interactions of gauginos 
f f M^ — 

/ d^2/£4+D=/ rf^y— Vtr(vl/zr'^K,M/]) ^ Ax0X, (107) 

JXd J Xd ^9* 

with the effective Yukawa coupling constants given by 



' Xd 
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Thus the low-energy Yukawa couphngs are given by overlap integrations, just as in fl22l) . where 
now is interpreted as the gauge coupling constant of the super Yang-Mills theory. 

Having specified the origin of the quarks, leptons, and their Yukawa couplings, we know what 
the underlying gauge group G has to be. For H = SU(5)gut5 super Yang-Mills interactions of 
G = E-i or Eg give rise to all of the Yukawa couplings of the quarks, charged leptons and neu- 
trinos [19]. For the Pati-Salam group H = SV{4)c x SU(2)l x SV{2)r, G = SU(8) is sufficient. 
When the Heterotic Eg x E'^ or SO (32) string theory is compactified on a real six-dimensional 
compact manifold, we obtain Yukawa couplings in this way; for example, with G = Eg and a 
six-dimensional manifold Xd=6 for the Eg x Eg string theory. In some compactifications of the 
Type IIA string theory (or supergravity on eleven dimensions) and Type IIB string theory (or 
F-theory), the Yukawa couplings of quarks and leptons arise from open string interactions. For 
example in the Type IIA string theory with an intersecting D6-D6 system, with quarks and lep- 
tons localized at the intersections of D6-branes, the Yukawa couplings are generated by spanning 
a world sheet to three intersection points with D6-branes connecting them. However, some of 
these compactifications can be understood as a certain limit of Heterotic compactification; this 
is the essence of string duality. Thus the toy models that we use to calculate Yukawa couplings 
can capture certain aspects (and maybe all) of these compactifications as well. Therefore, the 
numerical analyses of this paper can be used to understand the statistics of fiavor structure 
arising from a large region of the string theory landscape. 

If quarks, leptons and the Higgs originate from a super Yang-Mills theory of G containing H, 
and the Yukawa couplings from the super Yang-Mills interaction of G, then we cannot arbitrarily 
assume the following: 

• the shape of zero-mode wavefunctions — these must be solutions to the equations of motion. 

• the number of zero-modes — this is determined by the topology of the geometry of the 
internal manifold and of the gauge field configuration on X^. 

• the number of independent parameters that freely scan — only moduli parameters, i.e. de- 
formations of the manifold and gauge field configuration that cost no energy, are scanned. 

The possible restrictions above were ignored in the preceding sections, as we introduced assump- 
tions instead. In this section we discuss when such assumptions are justified, when they are not, 
and when not what one should expect instead. In sections 17.11 and 17.21 we provide a pedagogical 
and bottom-up introduction to ideas that motivated us to assume Gaussian zero-mode wave- 
functions with center coordinates that scan over the landscape. Most of the content of these 
sections is not particularly new, and the intended audience is non-string theorists. We use the 
remainder of this section to describe what the toy models of sections [3], [5] and [6] — i.e. Gaussian 
landscapes — mean in the context of the string theory landscape. 

Before we proceed, let us comment on the basis-independence of observables. The Standard 
Model has three independent fermions in a given representation, and the 3x3 Yukawa matrices 
transform under a group U(3)g x \j{3)u x U(3)j x 11(3); x U(3)e that preserves the canonical 
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kinetic terms of the fermions. Observables such as mass eigenvalues and mixing angles also do 
not change under this transformation. Thus, any probability distribution of Yukawa matrices in 
a landscape-based theory should be either invariant under the U(3) transformations, or defined 
only for classes of Yukawa matrices equivalent under the U(3) transformations. This basis 
independence follows in landscapes that derive from super Yang-Mills interactions. Indeed, 
the zero modes — i.e. the solutions of massless Dirac equations in a given representation — form a 
vector space, the rank of which gives the number of generations. One should find an orthonormal 
basis of the vector space {ipi}, such that the metric of the kinetic term is given by 



and use the zero-mode wavefunctions of the basis vectors to calculate the elements of Yukawa 
matrices through 



The above U(3) transformations correspond to basis transformations of the vector space of zero 
modes. Here the U(3) basis transformations correspond to no more than choosing different de- 
scriptions of the same vacuumH and the actual statistical elements of the landscape correspond 
to vacua, not Yukawa matrices. Thus, any landscape generating vacuum statistics through the 
overlap integrals of zero- mode wavefunctions is basis independently 

The Gaussian landscapes introduced in sections [3l [5] and [6] are meant to be simplified versions 
of (subsets of) the landscape that derives from string compactification. Because of the way we 
have simplified the landscape — scanning the center coordinates of zero modes completely ran- 
domly and independently — the basis independence of the landscape from string compactification 
is lost in our Gaussian landscapes. However if the Gaussian landscapes are regarded as tools to 
acquire a better understanding of flavor physics in the landscape from string compactification — 
and this is certainly our attitude — then the lack of basis independence in Gaussian landscapes 
is itself not a major problem. 

7.1 Domain Wall Fermion and Compactification of Field Theories 
on Six Dimensions 

Domain Wall Fermion 

^^Only the vector space composed of solutions to the zero-mode equations of motion have these U(3) "symme- 
tries;" the interaction terms — which are trihnear (not quadratic) in fields on a given background — do not respect 
these flavor symmetries at all. 

^^The U(3)-invariant Harr measure for mixing angles follows in basis independent landscapes [13] when A" 
and A'' (A*^, X'^, and Ca/s) are independent. In landscapes derived from overlap integration, however, there are 
correlations between the relevant matrices, and hence the mixing angles do not follow the Harr measure. 




(109) 




(110) 
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It is well known that a chiral fermion in a four-dimensional effective theory is localized at a 
domain wall in a five- dimensional theory [2U] . Consider a fermion in a five dimensions, 



d'xdx,C5 = d'xdx5 ix,^) [i{7''d, + ^'d,)-<P{x,)] ^ ^ , (111) 

where 'j^ = i on ip and —i on x- A chiral zero mode is localized about the point X5 = x^^ where 
a background scalar field configuration (/'(xs) equals zero, its wavefunction being 

When the background configuration is approximated by a linear function ^(xs) ~ F{x5 — x^^q), 
then the zero-mode wavefunction is approximately Gaussian: 

ij{y) ^ e-^("«-"^'°)'/^ (113) 

Thus the Gaussian profile follows from the zero-mode equations of motion rather easily!^ 

However, when the x^ direction is compactified into S^, the scalar field ^(xs) must be periodic. 
If ^(xs) hits zero with a positive slope k times along S^, it does with a negative slope k times 
as well. Thus, when there are k chiral left-handed zero modes in the low-energy effective theory, 
there must be k chiral right-handed zero modes as well, and the net chirality is zero. Furthermore, 
the background field configuration (f){x^) should satisfy its equation of motion, but kink-anti-kink 
configurations are in general not stable. 

Domain Wall Fermion from a D = 6 \J{1) Gauge Theory 

These problems are addressed if ^(xs) is not a scalar field. Consider for example a U(l) gauge 
theory on a six-dimensional spacetime with a charged fermion: 

£6 = ^zr^DM*, DM = dM-tAM, (114) 

where the gamma matrices are chosen to be 

rM=o,i,2,3 = T''=^'' = 7^ ® . (115) 

When the gamma matrices are chosen as above, the S0(5,l) spinor \1/ consists of four four- 
dimensional Weyl spinors, = [ip^jXp V^j.) Xj.)"^- ^ basis where r°^^^^^ = 1 ® and 

= 7^^ ® r\ = 7^ ® F^ = 1 ® ^r^ (116) 

the six-dimensional Dirac equation is given by 



^^On a domain- wall background, the zero-mode wavefunction is Gaussian in the middle of domain wall, where 
(xs) is approximately linear, and decreases exponentially outside the wall, where ^(xs) is nearly constant. 
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where D = {d^ — lA^) — i^d^ — iAq) and D = {d^ — iA^) + i^d^ — iAq). Let us focus on the 
Weyl spinor When is ignored, this spinor's Dirac equation is the same as the five- 

dimensional equation above — that is Aq enters the equation just as ^(xs) does in fillip . Thus Aq 
acts hke a mass parameter for a five dimensional fermion. Chiral fermions are localized where 
Aq "vanishes," just like with the domain wall fermion. Below we make more precise what we 
mean by "vanishing" Aq. 

Let us consider a U(l) gauge theory compactified on T^, with coordinates X5 and xq having 
periods L5 and Lq. Suppose that Aq is linear in and A5 is constant: 



^5 = ^5 



Aq = FX5 + Cq . 



:il8) 



The field strength F^q = F is constant, and satisfies an equation of motion dmF^^ = 0, where 
■m,n = 5, 6. There is no issue of kink-anti-kink instability (however see the discussion at the end 
of this section, where another instability is discussed). Since Aq is part of a gauge field, it must 
be periodic in X5 — + L5 only up to a gauge transformation. Since a gauge transformation 
ip{xQ) = e^mi's/ie shifts Aq by 2tx / Lq., Aq{x^) and ^6(^5 + L5) can differ by an integral multiple 
of 2ti / Lq. Therefore FL^ (and the field strength F) have to be quantized: 

2'Kk „ 2'Kk 



FU 



It follows from the index theorem that there are 



L^Lq 



kez. 



k= - 



F 

rp2 2TT 



(119) 



(120) 



chiral fermion zero modes. Because of the linear configuration of the effective mass term Aq, 
zero modes correspond only to left-handed spinors (if k > 0), and net chirality is obtained in 
the low-energy effective theory. 

Fermion Zero-Mode Wavefunctions: I 

It is easy to find one of the k zero-mode wavefunctions referred to above: 

i/j^{x5) ^ e^^^^^e"-^ (^^+^)' , Xi{x5) = , (121) 



just like a domain wall fermion. This zero mode is localized where the effective mass parameter 
Aq = Fx5 + ^6 is zero; x^^ = —^q/F. Since a gauge transformation (p{xq) = e27rifcx6/L6 shifts 
Aq by 2Tik/LQ and shifts —S,q/F by — the center coordinate x^^ of the Gaussian profile can 
always be chosen within the interval [0,^5]. To obtain the other {k — 1) zero modes, note that 
a Kaluza-Klein momentum in the Xq direction effectively shifts the Wilson line ^g- 



(^6 - e^"^^6 %1}{x^,xq) 



Oq 



^6- 



27rp 

17 



iP{x5,Xq) . 



(122) 
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Hence the remaining zero modes are 

^t(x5,X6) ^ e^«5"«e'"^^e-*("^+^^^^0' , Xi{x5,xe) = 0, (123) 
for p E 1i (mod k). The center coordinates of these k Gaussian zero modes are at 

X5 = + for p e Z/kZ. (124) 

The width of the Gaussian profile is 



d=^ = \R^- (125) 

It has been assumed that the width of Gaussian wavefunctions d can be parametrically smaller 
than the size of the extra dimension in the direction. This is equivalent to 

d 1 



-^5 \J'lT\ku ' ^ ^ 

where u = L^/Lq parametrizes the shape of T^. In the analysis of section [3] (the Gaussian 
landscape on S^), we found d/L^ ^ 0.08 was sufficient to reproduce the hierarchy of the quark 
sector. This translates to -u ^ 8 when k = 3. 

Fermion Zero-Mode Wavefunctions: II 

In fact the fermion zero-mode wavefunctions on are known exactly; the approximate form 
fll23p is valid only when can be ignored, i.e. when Lg is so small that only one Kaluza-Klein 
momentum is relevant at a time. Let us take a detour here to see the form of the zero-mode 
wavefunctions when L5 ^ Lq is not satisfied. First, we note that the fermion obeys a twisted 
periodicity condition. 



7/^1 (x5 L5,X6) = e ^V'T (2^5, are) , 'il:^{x5,xe + Lq) = ip^{x5,xe) , (127) 



just like. 



dm - iAm{x5 + L5, xe) = e^''*''^6 [d^ - iA^ix^, Xq)] e ^'"'''4 , (128) 

dm - iAm{x5, Xe + Lq) = dm- iAm{x5, Xq) . (129) 

The wavefunctions of the k chiral zero modes, which satisfy iDtp-^ = 0, are given by [21] 

^^{x,,xer = C^e-2-*fV«^-«e'^("^+^)("'^+^^'^-^)e-^("«+^^«^^)' (130) 

for p G Z/kZ = {0, 1, 2, ■ ■ ■ , (/c — 1)} (mod A;). Here we have defined 

= + — , = xe - — , v = , T=-r] To=i— = - . (132) 

F F Ls k L5 u 
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In addition, 'i9o;0 is a theta function defined by 

mez 

The wavefunctions (11311) form an ortlionormal basis wlien the normahzation constant is 

C=(^]' ^ (134) 
\kuj y/M^L^Le 

If Lq ^ L5 and Lq ^ 1/ \/F ^ \JL^L^^ then the zero-mode wavefunction (11301) receives dominant 
contributions only from the k consecutive values of m for which X6,o = F —{m/k)LQ is between 
and Lg. For other values of m the last factor is suppressed by a power of e"*^'^/^-'*-^''/'^'^^ <^ 1 
when xq G [0, Lg]. Thus the wavefunctions are linear combinations of Gaussian profiles localized 
in the direction, centered at xq^ = — {m/k)LQ. 

If L5 ^ Lq, then the expression given in (11301) is not useful to obtain an intuitive picture. 
Using a modular transformation of the theta function, 

^0;o(t'; r) = z^e-'^^'^o-^o (j; , (135) 

(I13ip can be rewritten 



xp{x5,xey = Ce'«^^^e-^(^5)'e-^(*5+^^«)' X iv^e~^(*^+**«-^^^) eo;o(t'';r') (136) 
= Czv^e^«^^^e'"'^^e-^(^^-^^^)'i9o;o(^^';r') , (137) 



kL^ f X5 + ixQ p\ , 1 k 

T = — = = iku . 138 



where 

^ Hq \ kj ' ' ~ T To 

When M ^ 1, all but one term in the series expansion of the theta function (I133P are suppressed 
by powers of e"'^'^" <^ 1, while the leading order term is a Gaussian wavefunction localized at 
X5,o = —i^/F + {j>/k)L^. This is the solution we obtained in (11231) . 

Although we have chosen a particular gauge (11181) . the localization pattern of fermion zero 
modes does not depend on the choice of gauge: the U(l) gauge transformation on fermion fields 
is not more than a phase multiplication. Although we chose Aq to be X5 dependent while 
is not xq dependent, the fermion zero modes exhibit a localization in x^ when Lq <^ L^, and a 
localization in xq when L5 -C Lq. Thus the width of the Gaussian profiles is given by (I125P in both 
cases, and is smaller than L5 and Lq, respectively. As expected, despite the x^-xq asymmetric 
gauge choice, the localization pattern is not x^-Xq asymmetric. The modular transformation 
property of theta function guarantees this. 

Note that when L5 and Lq are comparable, the zero-mode wavefunctions are neither Gaussian 
nor localized. For more about zero-mode wavefunctions on tori, see for example |22j . 



65 



Zero Modes of Various Fields in Different Representations 



At the beginning of this section, we advertised the idea that all of the fermions, the Higgs boson 
and the gauge bosons are unified into a super Yang-Mills multiplet of a gauge group G that 
contains the gauge group of the Standard Model or some unified theory. Let us describe this 
idea a little more explicitly, using the language of compactification with a U(l) gauge field 
background. Consider a G = U(6) super Yang-Mills theory (a parallel description applies to 
Pati-Salam unified theories with G = SU(8) and H = SU(4)c x SU(2)l x SU(2)/j). Among the 
generators of the U(6) symmetry, U(l) field strengths on are introduced along 

ts = diag(l,l,l,0,0,0) , = diag(0, 0, 0, 1, 1, 0) , ti = diag(0, 0, 0, 0, 0, 1) . (139) 

SU(3) X SU(2) and a couple of U(l) factors commute with these field strengths, and these are 
identified with the Standard Model gauge group. Gauginos in the 3x2 block and in the 1x3 
block are in the (3, 2) and (3, 1) representations of SU(3)c x SU(2)l, just like the quark doublets 
and the anti up-quarks. Thus, the corresponding zero modes of the gaugino can be identified 
with these Standard Model fields. Meanwhile, the Higgs boson may arise from Am in the 2x1 
block. Then the gauge interaction of the gauginos in fill 01) gives rise to the up-type Yukawa 
coupling: following the commutation relations of U(6) (c.f [23]) we have. 
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uqHu. 



(140) 



The Yukawa couplings are then calculated by overlap integrals of the zero modes. 

The U(l) fiux in U(6) is given by tgF^^) + taF^^) + tiF^^), with each of F^") = ^5^^"^ 
(a = 1,2,3) given by A'^^^ = (27rA;("V^5^6)a;5 + as in ffTT8|) and ffTTOD . One has to set the 
field strengths F^^\ F^"^^ and F^^^ so that there are three quark doublets and three anti-up-quarks 
in the low-energy spectrum. This means that 



ir{3) ^(2) 

rp2 1T\ 1T\ 



2^2 27r Itx 



A;(^) - = 3. 



(141) 



The three Gaussian zero-mode wavefunctions of quark doublets and anti-up-type quarks are 
localized at 



X 



+ 



3 



27r(A;(i) - A;(3)) ' 



(3)^ 



+ 



J G Z mod (A;(^) - A;^^)) , (142) 
i e Z mod {k^^^ - k^^^) . (143) 



This is an illustration of the picture described in the introduction to this section. We started 
from U(6) super Yang-Mills theory on and broke the U(6) symmetry by turning on a gauge 



66 



field background, so that an SU(3)c x 811(2)^ gauge symmetry is left unbroken along with some 
U(l) factors. Furthermore, the gauge- field background generated chirality in the low-energy 
spectrum, and determined the shape of the zero-mode wavefunctions (which are Gaussian when 
L5 ^ Lq or Lg ^ L5). Wilson lines ^^""^ (a = 1, 2, 3) are chosen arbitrarily, yet the field strength 
satisfies the equation of motion. Thus, these constants determine the zero-mode wavefunctions, 
and hence the observables that arise from the Yukawa matrices. If the Wilson lines are scanned, 
the peaks of the Gaussian zero modes are also scanned, and hence the Yukawa couplings are 
scanned. Since the zero modes of quark doublets and anti-up-type quarks depend differently 
on the Wilson lines — see fll42p and (11431) — the peak positions of the wavefunctions of different 
representations are scanned independently from one another. 

Problems 

The compactification of super Yang-Mills theory reproduces certain aspects of the Gaussian 
landscape on S^, but it also has serious problems. The easiest problem to identify involves the 
sum rule on the chirality. Because of the simple relation {F^^^ — F^"^^) + {F^^^ — F^^^) + (F^^^ — 
F^^^) = 0, we have 

i^k = -(#g + #m) = k^'^ - k^'^ = -6 , (144) 

implying that the low-energy spectrum contains n > up-type Higgsinos and 6 + n fermions in 
a representation conjugate to that of the up-type Higgsinos. One encounters a similar sum rule 
in any compactification on a smooth two-dimensional manifold with U(l) field strengths; thus 
this is too simple a framework to provide both a realistic spectrum and Yukawa interactions. 

A more serious problem is that the constant U(l) x U(l) x U(l) gauge field configuration 
in U(6) is not stable. Although this configuration satisfies the equations of motion, this only 
means that it is a stationary point of the action. In fact, the equations of motion of the gauge 
field in the 3x2 and 1x3 blocks are given by 

[A,, + {F('^ - F^%e] {A, + ^A,) = 0, (145) 

[A56 + (F(l^-F(3))56](A5+2A6) = 0. (146) 

Here A56 is the Laplacian on the x^-Xq plane. The positive (F^^^ — F^'^^)^q and {F^^^ — F^^^)^q 
required for the proper quark chirality implies that the low-energy complex scalar fields coming 
from A^+iAg have negative mass-squared. Therefore the constant U(l) gauge field configuration 
is not stable on T^. It is stable only when F^f^ = F^ = F^l\ which yields no net chirality. 

Complex Valued Wavefunctions 

In section 17.21 we move on to discuss more complicated compactifications of supersymmetric 
Yang-Mills theory in search of frameworks for constructing Gaussian landscapes that avoid the 
above problems. For the remainder of this section, however, we digress to explore how CP- 
violating phases can be introduced into these landscapes. We look at the simplest geometry 
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imaginable — compactification on — considering that the essence of obtaining CP-violating 
phases will not be lost in more complicated frameworks. 

Real-valued Gaussian wavefunctions never yield complex phases in the Yukawa matrices or 
CP-violating phases in the quark and lepton mixing matrices. Although the zero-mode wave- 
functions have phase factors Q^^^^^Q'^'^^vi^a/La) ^ these do not end up as complex phases in the 
Yukawa matrices. The first factor e*^^"^^ cancels in the integrand of the overlap integration. For 
example, consider the U(6) SU(3)c x 811(2)^ x U(l)y symmetry breaking model. The up-type 
Yukawa operator is neutral under any one of U(l)^"-' (a = 1,2,3) symmetries, and this ensures 
that the e**5 phase factor in the quark-doublet wavefunctions is canceled by the factor e ^'^ 
in the anti-up-quark wavefunctions. As for the second phase factor Q'^T^wi^e/Le) ^ recall that the 
overlap integration is carried out on both the X5 and Xg coordinates. This phase factor is a 
plane wave associated with Kaluza-Klein momenta, and these momenta are conserved in 
compactifications. Therefore in any matrix element give n by overlap integration on T^, these 
plane wave phase factors cannot yield complex phaseso 

However, there is a simple way to obtain complex phases in Yukawa matrices without con- 
sidering more complicated geometries. So far we have assumed is rectangular; that is the two 
independent periods, 

(X5,X6) ~ (X5 + L5,X6) , (xg, Xe) ~ (xg, Xg + Lg) , (147) 

are rectangular. In other words, using the complex coordinate v = (xs + ixQ)/L^ we have, 

V ^ V + 1 , V V + To , (148) 

with a pure imaginary tq given by (11321) . When tq is not pure imaginary (while keeping the 
metric ds"^ = dx\ + dx\), zero-mode wavefunctions have more complicated complex phases, as 
we see explicitly below. 

We now use the gauge field configuration, 

^5 = ^5 , A^ = F ( X5 - -X6 ) + ^6 + -^5 = f ^5 - -^6 ) , (149) 

where Ti and T2 are the real and imaginary parts of t = ro/k, and 

- (150) 



Lpm(ro) LiT2 



^^Kaluza-Klein momenta conservation leads to an approximate texture in the Yukawa matrices. When only the 
leading term in the series expansion of the theta function is kept for each of q, u, and /i, a given quark doublet has 
a non-vanishing Yukawa matrix element with only one anti-up-type quark. All of the other matrix elements are 
suppressed by at least e~'^", which arises after sub-leading terms in the series expansion of the theta function are 
included. On the other hand, when compactifying on generic manifolds there is not necessarily any Kaluza-Klein 
massless U(l) gauge field and its corresponding conserved Kaluza-Klein momenta. This is why we consider that 
the Yukawa texture from Kaluza-Klein momenta conservation is an artifact of torus compactification, and have 
ignored it in the Gaussian landscapes of earlier sections. 
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The fermion fields satisfy the twisted periodicity conditions in (11271) for the periods in (I148p . 
The zero-mode wavefunctions (I13ip are modified and become 

= e^«=^^e-^(^)'ev^^^O;0 {v - |; r) , (152) 

where v and f are the complex conjugates of v and r, respectively. To see the behavior of these 
wavefunctions when T2 = Im(ro) is small (roughly equivalent to m ^ 1), we use fll35p . Dropping 
the sub-leading terms in the series expansion of the theta function and maintaining only the 
0{t2) and 0{1) terms in the exponent, we have 



<(X5,X6) ^ e*««^'^e"^" Vxt-i; f^e^ . (153) 



The factors e*^^^^ and e^'^^'^a*-^/'^-**-^^/^^'' may cancel in overlap integration. The second factorQ 



^■2 



has a complex phase when r = —t[/t2 7^ 0. This phase will not cancel with overlap integration 
on T^, which is why we used this second factor in the zero-mode wavefunctions in the Gaussian 
landscape of section 16. 4[ 



7.2 Doubling — x Compactification of Field Theory on Eight 
Dimensions 

We have seen that the zero-mode wavefunctions are localized and Gaussian in a compactifica- 
tion with a constant U(l) gauge field background, when L5 ^ Lg or vice versa. Thus, when the 
is doubled and = x compactification with a constant U(l) background is considered, 
the zero-mode wavefunctions are localized and Gaussian in two of the four directions, when 
is short in the other two directions. As we will see in this section, in compactification the net 
chirality of various representations is no longer subject to a linear sum rule, and the instability 
problem in the compactification can be avoided when the gauge field configuration satisfies 
an anti-self-dual condition. Note that sections 17.21 and 17.31 discuss supersymmetric Yang-Mills 
theory on eight dimensions as a possible origin for Gaussian landscapes, but not because we 
regard a theory on eight dimensions as a candidate for a fundamental theory. Instead, these 
sections are intended to provide a step-by-step introduction to ideas that formal theorists are 
already familiar with, but with minimum technicality. Later, in sections 17.41 and [73| we consider 
a realistic framework for building Gaussian landscapes. 

Fermion Zero Modes in Anti-self-dual Gauge-Field Background 



30Here we have used 27rT^ ~ {L^/df; see and 
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We label the four coordinates of x^^j^s and the circumferences in these directions are denoted 
-^5,6,7,8- We assume that these directions are all orthogonal, i.e. the internal metric is 



ds^ = {dx^Y + {dxfsf + {dxjf + (c/xg)^ . (155) 
On this x T^, we introduce a U(l)2 gauge field configuration embedded in SU(2): 

I (2) _ f{2)^3 4(2) _ / 27rA; (2)\ 3 



At' = e,'r^ , At' = [j^x, + it^j , (156) 

4^^ = -efV, 4^^ = -(|^-r + e^)r^ (157) 

The field strength is constant and satisfies the equation of motion. Moreover, this configuration 
is stable when the ant i- self- dual condition F = — * F is satisfied; that is when 



2tt k 271 k' 



56 — —-^78 ! 



'158) 



A Weyl fermion in the fundamental representation of the SU(2) has 2kk' zero modes on T^. 
To see this, we use the following gamma matrices: 

rM=o,i,2,3 _ ^ 1 ^ 1 ^ = 7^ ® ^ 1 , = 75 ^ ® . (159) 

A Weyl fermion on eight dimensions consists of four- dimensional Weyl fermions (V^fj, V^Jp xft; Xu)) 
where a is a "color" index running a = 1, 2 in the case of a fermion in the fundamental repre- 
sentation of SU(2). Corresponding to the kk' different choices of p G Z/A;Z and p' G Z/A;'Z, the 
Dirac equation on with the gauge field background (11571) has kk' zero modes in ijj"-^^: 

^}j1='-'P'P' = ^P^(x,,x,;kr^|J^ixr,Xs;kr\ ^ J^^^'^' = , ^'^= W = q , x = 0, (160) 

Likewise, there are also kk' zero modes in tp"-^"^: 

riT'''"'"' = M^5,x,;kr^^{xj,Xs;k'r' , C'^'''' = 0, r=''''''' = 0, x = 0. (161) 

When L5 ^ Lq and Lj ^ Lg, all of these wavefunctions are approximately Gaussian in the 
X5-X7 plane, since the wavefunction is a product of Gaussians in the X5 and X7 directions. Two 
zero modes, one in ipfj"^ and the other in ipf^^, are localized at each of 



The width of all of these zero modes is d = 1/ a/ | -F56 1 in the X5 direction and d' = 1/ ^/\Fts\ 
in the Xj direction. The anti-self-dual condition (11581) implies that the two widths d and d' are 
equal, as was assumed in the D = 2 Gaussian landscapes of section [5l 

Absence of a Chirality Sum Rule 
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The two serious problems of compactification with a U(l) gauge field background derived 
from instability of the gauge field configuration and the sum rule satisfied by the chirality of 
different representations. The anti-self-dual condition fll58p addresses the first issue in (or 
in any four-dimensional manifold) compactification, while in what follows we see that there is 
hope to resolve the second issue in four-fold compactifications as well. 

The chirality sum rule on derived from the Yukawa interactions being consistent with 
the gauge symmetry of an underlying group G. Therefore, to discuss the existence of such a 
sum rulqfj on T^, we must consider constraints from the underlying gauge symmetry whose 
super Yang-Mills interactions give rise to the Yukawa couphngs of quarks and leptons. Here we 
consider G = Eq and H = SU(5)gut- Eq contains a subgroup 

Ee D SU(2) X SU(6) D SU(2) x U(l)6 x SU(5)gut • (163) 

The H = SU(5)gut symmetry is unbroken when a background gauge field on is contained 
within SU(2) x U(l)6- We use the U(l)2 gauge-field configuration Am in (11571) . embedded in 
the SU(2) factor, and introduce a U(l) gauge field background in the U(l)6 factor, given by 

(6) _ 2nk" (6) (6) _ f 2TTk"' (g) 



= X5 + , = - "j^xr + , (164) 

with another anti-self- dual condition 2Txk!' jL^L^ = 2T[k!" / L-jL^. 
The irreducible decomposition of the Eq Lie algebra, 

ee-adj. (adj., 1) + (1, adj.) + (A^G, 2) , (165) 

(adj.,l)° + (5,l)^ + h.c. + (l,adj.) + (A'5,2)-3 + h.c., (166) 

shows that g/f) contains both (A^5, 2)"'^ and (5, 1)^, candidates for the 10 = (g, m, e) and H{h) 
multiplets, respectively. The above decomposition refers to the the SU(6) x SU(2) subgroup in 
the first line, and SU(5)gut x U(1)6 x SU(2) in the second line. Gauge interactions of eg gauginos 



■^^To be precise, chirality on four dimensions is noi obtained from a super Yang-Mills theory on eight dimensions 
compactified on with an anti-self-dual gauge field configuration. This is because the super Yang-Mills multiplet 
contains a pair of Weyl fermions with opposite chirality. When a zero mode is found in V'^j^^ (^if^) ^ Weyl 
fermion ("011: V'i.T: XTT' Xii)j another Weyl fermion with the opposite chirality, (V'Ti: V'iTi XTT' Xii)j also has a zero 
mode ^5^^ (V^J^^). The chirality we refer to in this section concerns compactification of ten to six dimensions, 
which has nothing to do with the chirality on four dimensions. As we explain later, however, the chirality on 
six dimensions in compactification has non-linearity, which is shared by chirality on four dimensions in any 
six-dimcnsional compactification. This "chirality" is discussed in the context of compactification because this 
is the simplest system with this non-linearity. 
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may in principle give rise to the four- dimensional up-type Yukawa couplings 

C4 = ipio ■ ipw ■ <pH{5) ■ (167) 

As a quick check, note that in the ee algebra tr ((A^5, 2)"'^ [(5, 1)"^^, (A^5, 2)"'^]) does not vanish. 
The coefficients (Yukawa matrix elements) are calculated by the overlap integration of the zero- 
mode wavef unctions, picking up a pair of V^^]"^ and V'j'f^- If compactification preserves 
A/" = 1 supersymmetry, a complex scalar and a chiral fermion in a chiral multiplet have the 
same wavefunctions on the internal manifold. Therefore the wavefunction of the Higgs boson is 
Gaussian, whenever the Higgsino wavefunction is Gaussian. 

The number of zero modes is calculated in terms of the flux quanta k, k', k" and k'" . Up-type 
Higgsinos in low-energy spectrum are zero modes of the Dirac equation involving a U(l) gauge 
field GAm''. Therefore the number of up- type Higgsino zero modes is 

= {6k"){6k"') = 36k" k'" . (168) 

Light fermions in the 10 representation are zero modes of a Dirac equation with gauge field 
A^'^^ — 3A^^^ . The number of zero modes is 

#10 ={k- 3k"){k' - 3k'") + {k + 3k"){k' + 3k'") = 2kk' + 18k"k'" , (169) 

with the ffist term coming from those in iIj"-^^ and the second term from those in ip"'^'^. We have 
four (discrete) parameters to choose by hand to fit the desired multiplicities of the fields in the 
two different representations, and 10. 

The discussion so far has not been realistic. The above multiplicities count only the num- 
ber of hypermultiplets of A/" = 2 supersymmetry in four dimensions. We have also commented 
in footnote [3S1 that the Higgs scalar field in the up-type Yukawa couplings should belong to an 
M = 2 vector multiplet. However, it will now be easy to take one more step and triple to a 



■^^Let us first consider how gauge indices are contracted. The SU(5)gut indices of the two fermions V-To ^re 
contracted with the SU(5)gut totaUy anti-symmetric tensor, and the contraction is symmetric under the exchange 
of the two. On the other hand, the internal gauge indices of SU(2) are anti-symmetric as the above interaction 
is neutral under the underlying gauge symmetry SU(2) x U(l)6: when the contraction is in the first term of the 
right-hand side oi2~^ ^2~^ ^ 1+^ ~ A^2 • • • . Generally, gauge indices are contracted anti-symmetrically under 
the exchange of the two zero modes in the (A^5, 2)"'^ representation. This should be the case because the gauge 
indices are contracted through the structure constant of a Lie algebra in piOp . 

Spinor indices are contracted symmetrically under the exchange of V'lo's; both the SO (3,1) contraction of 
left-handed spinors in four dimensions and the contraction of spinor indices of SU(2) C S0(4) in the internal 
space are ant i- symmetric. Therefore the combined contraction of gauge (anti-symmetric) and spinor (symmetric) 
indices is anti-symmetric. This is expected since the gauginos in eight dimensions are Grassmann variables. 

The Higgs scalar should, then, originate from the complex scalar of a super Yang-Mills multiplet on eight- 
dimensions. However, because of the A/" = 2 supersymmetry, preserved in compactification with an anti-self- 
dual gauge field configuration, the existence of massless Higgs boson implies that its A/" = 2 super partner, a vector 
field in the same representation, should also be massless. This implies that the symmetry is not broken down to 
SU(5)gut- As we have seen, T'^ compactification with an anti-self-dual configuration has so many problems that 
it cannot be regarded as a realistic framework for landscapes. Our goal in section [7721 is to introduce field-theory 
ideas that apply to many compactifications, using the language of compactification. 
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compactification of a super Yang-Mills theory in ten dimensions. Chirality on four dimensions 
is obtained in D = 6 compactifications, and the net chirality of a given representation is cubic 
in discrete parameters specifying the fluxes on T^. No sum rule holds among the chirality of 
various representations, since they are not linear functions of the flux parameters. As more flux 
parameters are involved, it is easier to fit the multiplicities of fields in various representations. 
It is also known that one of the gauge background stability conditions is F^^g"^ = (in com- 
pactifications preserving A/" = 1 supersymmetry) , and the compactification has instability 
because only one term contributes to the left-hand side, with no chance of cancellation. Hence 
this instability is an artifact of the D = 2 compactification. Eg contains E^, and therefore the 
algebra that led to the up- type Yukawa couplings still works in E^. Thus super Yang-Mills 
theories of ten dimensions compactified on six-dimensional internal manifolds can be a good 
theoretical framework for building Gaussian landscapes. 

Scanning of the Center Coordinates 

Before closing this section, let us see how the localized wavefunctions behave as the constant 
gauge field background on changes. We already know that all the zero-mode wavefunctions are 
approximately Gaussian if L5 ^ Lq and L-j ^ Lg. The center coordinates of these wavefunctions 
depend on the Wilson lines ^g^'*'*'^'* and C,^^'^^\ according to 

(x^,x^) = ff^-i^jL,,f|p^-^|L. 



I X5 , Xj 




10 ION _// Pit (C-3ef)L6\ / pTo' (d'^-3ef^)L 



)a=l h — WI O-rrih — '\hii\ \ ^' 



2Ti{k-'ik") J "'U'-SA;'" 2TT{k'-3k" 

{X^ ,Xj )a=2 




2Ti{k + 3A;") / \k' + 3k'" 2TT{k' + 3k'") 



where the coordinates in the first line are those of up- type Higgsino(s), while the second and 
third lines are those of the 10 fermions, with the second and third lines coming from ip'^^^ 
and ip"'^'^, respectively. Note that the coordinates of the second and third lines show different 
dependences on the Wilson lines. Thus as the Wilson lines are scanned, the localized fermions 
in the 10 representation of SU(5)gut change their relative positions. We see that the scanning 
of Wilson lines allows for more than the center-of-mass scanning of Gaussian peaks to zero- 
mode wavefunctions in compactification. In the Gaussian landscapes we scanned the peak 
positions of Gaussian wavefunctions without considering any correlations or constraints among 
the various peak positions. The above situation is still far from this treatment, but so far we 
have only scanned the Wilson lines. In section 17.31 we see that there are other stable gauge 
field configurations, and as more gauge field configurations are scanned, correlation will be lost 
among the Gaussian peak positions. This is the rationale behind the absence of correlations in 
the Gaussian landscapes. 
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7.3 Instanton Moduli and Random Peak-Position Scanning 
Instanton Moduli 

We have now seen an explicit example of a stable gauge field background on parameterized 
by Wilson lines. In general, stable gauge field configurations allow for continuous deformations, 
such as Wilson lines, and the parameters of such deformations are called moduli. As for gauge 
field configurations on a four- dimensional manifold, we know that those satisfying the anti- 
self-dual condition satisfy both the Yang-Mills theory equations of motion and the stability 
condition. Such gauge field configurations, called instantons, are known to have numerous 
moduli parameters; the Wilson lines on are just a subset of the instanton moduli parameters. 
Therefore we now consider what happens when the instanton moduli parameters are scanned. 
Our primary interest is to explore the moduli space of stable gauge field configurations on a six- 
dimensional manifold, as this is relevant to the compactification of the Heterotic string theory. 
However, to warm up we first study instanton moduli on a four-dimensional manifold. 

The Eq — > SU(5)gut symmetry breaking model in section UTH generates the up-type Yukawa 
couplings. As long as a background configuration of Zq gauge field is contained in its su(2)-|-u(l)6 
subalgebra, the SU(5)gut symmetry remains unbroken. Note that the gauge field configuration 
does not have to be pure Abelian over the entire T'^ as is assumed in (11571) . The gauge field 
configuration A^^^ in (I157p can be replaced by any one of the SU(2) instanton configuration^ 
on T^. The 't Hooft solution is an SU(2) instanton configuration on a fiat Euclidean four- 
dimensional space [21]; it has 5/ moduli parameters, where the instanton number / is topological 
and defined by 



/ch2,fund.(F) = -|trf.„d. =-^/^^^'^^^)' 



(174) 



here fund, stands for fundamental representation, R is an arbitrary representation, and "ch2;_R" 
stands for the second Chern character in the representation R. The gauge-field configuration 
A^'^^ in (I157p has / = 2kk' instantons. Meanwhile, the 't Hooft solution. 



— ~Vmn'^n In 1 



(y - VjY 



(175) 



•^^ The number of fermion zero modes is given by the topology of the gauge field configuration and of the 
geometry, as in the case of compactification. It follows from the index theorem that the number of fermion 



zero modes in a representation i? of a gauge-field background is 



chRiF)A{TX) = - / ch2;fl(F) 

X J X 



dimi? 
24 



Pi{TX), 



(173) 



X 



where A{TX) is the A classes of X and pi{TX) is in the first Pontrjagin class. This expression is valid for 
an arbitrary four-dimensional manifold X. The first term on the right-hand side is equal to 2Ti^I, that is it is 
proportional to the instanton number /. As long as the instanton number remains the same, the net chirality 
does not change for any instanton configuration. See the text for the definitions of ch2, r{F) and /. 
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describes / isolated instantons, centered at {yj)m with a size pj. Here, a = 1, 2, 3 label the three 
generators of su(2) and is the eta symbol of 't Hooft. Thus there are indeed (4 + 1)/ moduli 
parameters. In addition to these five moduli parameters per instanton, there are three more 
moduli parameters describing how an instanton solution is embedded within the group SU(2). 
Thus, there are (roughly) 81 moduli parameters for the SU(2) J-instanton configuration. 

On a compact four-dimensional manifold X, the number of instanton moduli parameters of 
a gauge group G is given by [25] 

- 2 / ch2;adj.(^) td{TX) = ATgI - / C2(TX) + c^iTXf , (176) 

Jx Jx 

where Tq denotes the Dynkin index in the adjoint representation (also known as the dual Coxeter 
number), td{TX) is the Todd classes, and dimG is the dimension of G. When G is SU(A^), 
Tg = N and dim G = N'^ — 1. Apart from the second term, which is associated with compactness 
of the four-fold X, the first term reproduces the result of G = SU(2) — 8/ moduli parameters. 
The second term corresponds to an obstruction for lifting isolated / instantons into an anti-self- 
dual configuration over all of the X (c.f. |2S])- The Wilson lines 7 § on the compactification 
of section rr2] are part of the instanton moduli. However unlike the Wilson lines of T^, which are 
associated with a non-trivial 7ri(T^), most of instanton moduli here are not strongly associated 
with the specific geometry of T^. Thus the existence of instanton moduli is very robust. 

Wilson lines as Instanton Moduli Parameters 

So far it has not been clear how the Wilson lines on X = are related to the instanton moduli 
parameters, such as {yj)m and pm, in the 't Hooft solution fll75p : in particular the constant-field 
configuration does not have a well-defined center. To see explicitly that the Wilson lines actually 
correspond to instanton center coordinates, ADHM data [27| and the Nahm transformation [2H] 
are useful. 

The ADHM formalism [27] allows one to parametrize the moduli space of an instanton 
configuration on a fiat space M'^. A set of parameters called ADHM data describes all possible 
instanton gauge-field configurations and, in turn, the data can be extracted starting from the 
gauge field configuration. For simplicity, we consider only the /-instanton configurations of an 
SU(A^) gauge group. Among the data is an / x J matrix-valued A^. This part of the data is 
extracted by 

(ArnY' = {p'\xM = I d'xr'^x^r, (177) 

^ ^ Jr* 
where {p,p' = 1, ■ ■ ■ ? -^) are the zero- mode wavefunctions of a fermion in the fundamental 
representation of the SU(A^) gauge group. Since the zero- mode wavefunctions are determined 
by the gauge field configuration used in the Dirac equation on M^, the data Am carry some of 
the information of the gauge field configuration. It is known that the data corresponding to the 
't Hooft solution are 

iAm) ={yp)mSp'p. (178) 
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Thus the data Am provide a way to extract "center coordinates" from the instanton gauge field 
configuration. When one considers a family of instanton gauge field configurations by modifying 

Am{x) A'^ix) = Am{x) + U , (179) 

the fermion zero-mode wavefunctions change according to 

ipix) — > ^p'{x) = e'^^''"'i){x) . (180) 

Therefore the data {p'\xm\p) are also extracted by 

{KX'' = {pW\p) = {p'\ - idu\p) = [ d'xij'^'^d^J^ . (181) 

^ ^ Jr* 

The gauge field configuration (11571) is already in the form fll79p . and the zero-mode wave- 
functions of the fermions in the fundamental representation are also already provided. The 
corresponding data, now calculated by integration on T^, not M^, are ([21] and references therein) 

(iW.)"- (-£.o.-g.o)v.^ (-(€r)^»-«--(€F)^"°)v- 

(182) 

Therefore the Wilson lines y and Q can be regarded as center-of-mass modes of the instanton 
center coordinates in the appropriate limit of the instanton moduli space. 

Scanning over Instanton Moduli Space 

If the moduli parameters of the instanton configurations are scanned randomly (according to 
some measure), the gauge field configuration is determined for each choice of moduli parameter, 
as are the zero-mode wavefunctions. In the Eq — > SU(5)gut symmetry breaking model of 
section [7^ fermions in the 10 representation of SU(5)gut in the effective theory are zero modes 
of fermions in the representation of an SU(2) x U(l)6 instanton background. As the moduli of 
an SU(2) instanton are scanned, zero-mode wavefunctions of the fields in 10 = (g, u, e) vary. The 
center coordinates of these wavefunctions will be scanned almost randomly, since the instanton 
center coordinates can be chosen arbitrarily in (11751) . The zero modes of a fermion in the 
representation are not classified into the zero modes in and those in ip"-^"^. The instanton 
configuration is contained in a Cartan U(l)2 subgroup for only limited points in the instanton 
moduli space; for generic points in the moduli space the SU(2) symmetry is completely broken 
and the distinction between a = 1 and a = 2 is lost. In general each zero mode has non-zero 
wavefunctions in both a = 1 and a = 2. This means that there is no selection rule in the Yukawa 
couplings; if some zero modes were exclusively in a = 1 and all others were in a = 2, then the 
Yukawa couplings would have involved only a pair of zero modes, one from each group. 

The Gaussian landscapes of sections [31 and E] extracted these features. In particular, we 
applied these features to all of the fields in the Standard Model. Certainly the Higgs wavefunction 
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is determined only by a gauge field background of an Abelian symmetry U(l)6 in the — *■ 
SU(5)gut symmetry breaking model; therefore the complexity (and variety) of instanton field 
configurations of non-Abelian symmetries has nothing to do with the scanning of the Higgs 
field's center coordinate. But this is just an artifact of choosing G to be minimal for the up-type 
Yukawa couplings. For a larger underlying gauge symmetry G, such as G = Eg, the Higgs 
is regarded as a zero mode of a field in the A^5 representation in Es SU(5)gut symmetry 
breaking due to an SU(5) instanton. Gaussian landscapes are based on an expectation that the 
variety of stable gauge-field configurations is so rich that the random scanning of gauge field 
moduli results in (approximately) random and independent scanning of the center coordinates 
of the zero-mode wavef unctions. 

One will notice here that the three Gaussian zero-mode wavefunctions in a given represen- 
tation do not necessarily satisfy the orthonormal condition (11091) . when the center coordinates 
are chosen completely randomly. Our approach has been to try to implement the rich scan- 
ning of gauge field moduli in Gaussian landscapes at the cost of giving up basis independence. 
Since the orthonormal condition is violated especially when two center coordinates coincide, the 
probability distribution functions of the smaller eigenvalues of a Yukawa matrix in Gaussian 
landscapes may become unreliable as the approach the largest possible value, 2;^ ~ in the 
notation of this paper. Although the center coordinates of two instantons can coincide, the 
two zero-mode wavefunctions associated with them should be properly modified so that they 
remain orthogonal. This modification is not taken into account in the Gaussian landscape. It 
is a yet-to-be tested question in string theory whether the gauge field moduli are such that the 
random scanning of center coordinates is a relatively good approximation or not, and if not how 
the correct distribution functions of observables would deviate from the predictions of Gaussian 
landscapes. 

Non-Gaussian Wavefunctions 

We here note that zero-mode wavefunctions are not always Gaussian for arbitrary choice of 
gauge field moduli parameters. In fact, when the sizes of instantons pj in the 't Hooft solution 
(I175p are much smaller than the typical distance between the instanton centers \yk — yi\, the 
't Hooft solution is a collection of isolated BPST instantons. Fermion zero modes are localized 
around the instanton centers y ^ yj, and their wavefunctions decay in as a power of the distance 
\y — yj\, not exponentially. Unless the size parameters pj are extremely small, however, overlap 
integration using these non-exponential wavefunctions tends to be larger than the 10~^-10~^ 
required to match the quark and lepton Yukawa couplings of the first generation. 

There are also situations where zero-mode wavefunctions decay linear-exponentially, rather 
than as a power law. In section [7711 (17.21) . we discuss only (T^) compactification, where no 
topological 1-cycles (2-cycles) can shrink while keeping the volume of (T^) finite. However 
in more complicated geometries there are topological cycles that can do this. When a U(l) 
flux is introduced on such a topological 2-cycle, and if for some reason the 2-cycle shrinks. 
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then symmetry breaking by the U(l) flux can be locahzed in extra dimensions, c.f. [29]. This 
contrasts with the situation in sections 17.11 and \7.'2\ where the symmetry-breaking U(l) field 
strength is spread out homogeneously over the extra dimensions. If the symmetry breaking 
U(l) fiux is localized at a point in the extra dimensions, then fermion zero- mode wavefunctions 
decay linear-exponentially, with the exponent proportional to the distance from the symmetry- 
breaking source. Depending on the choice of parameters, the linear-exponential wavefunctions 
may or may not lead to as large a hierarchy as that which results from Gaussian wavefunctions. 
It would be interesting to study the effects on Yukawa-related observables of such a localized 
symmetry breaking source, but this is beyond the scope of this paper. 

To date, there has not been much investigation into what part of the moduli space of gauge 
field backgrounds is more statistically weighted in fiux compactification. It might be discovered 
that the statistic distribution is more weighted in regions that lead to Gaussian (or possibly 
linear-exponential) wavefunctions, so that large hierarchy among Yukawa couplings follows as 
a likely consequence. The other possibility is that environmental selection in favor of a light 
charged lepton (and/or quarks) enhances the statistical weight for such regions. We do not know 
which, if either, possibility is correct. However, for the landscape to account for the observed 
hierarchical patterns in the mass eigenvalues and mixing angles, we must assume one of these 
possibilities. As long as either of these possibilities is correct, then Gaussian landscapes should 
not be too terrible an approximation (or an effective description) of the landscape formulated 
by a super Yang-Mills theory on higher- dimensional spacetime. 

7.4 r^-Fibered Compactification and M (II A) -theory Dual 

In this paper we present numerical simulations of Gaussian landscapes involving only D < 2 extra 
dimensions. This is due to limited computational resources, not because there is any theoretical 
motivation to study landscapes involving D < 2. Indeed, more interesting would be to study 
Gaussian landscapes on D = 3, since such a model would directly simulate (some fraction of) the 
landscape of Yukawa couplings in string theory. A Gaussian landscape on a three-dimensional 
manifold B corresponds to Heterotic string theory compactified on a six-dimensional manifold 
that is a T^-fibration over B. To clarify, when one says that a m-dimensional manifold X is a 
T^-fibration over a manifold B, this means that the manifold X locally resembles T" x 
almost everywhere in X. For example, T"^ = x and T^(a;5^6,7,8) = ^^(^^e.s) ^ ^^(2^5,7) in 
sections 17.11 and 17.21 are trivial examples of = and T^-fibered geometry. Here B is called 
a base manifold, and T" a fiber. 

Suppose that a six- dimensional manifold X is a T^-fibration over a three-dimensional man- 
ifold B. We choose a coordinate system locally, so that the direction is parametrized by 
(xe, xs, xio) and a local patch of B by (xs, X7, xg). When the periods of in its three directions, 
Lq, Lg, and Lio, are small compared with the size of B (that is, the cubic root of the volume 
of B), then a stable gauge field configuration on X is approximately described by fields ^6,8, 10 
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that vary slowly on the coordinates x^^i^i^ of B. Fermion zero modes are localized at points on B 
where all of Aq{x^,Xt,x<^), A^{x^,XT,Xg) and y4io(x5, 0:7, Xg) vanish. As long as Lg, Lg, and Liq 
are all small, the zero- mode wavefunctions are approximately Gaussian [SU]. Although X is not 
globally the same as x or x T^, the local structure of X — its T^-fibration — is sufficient 
to determine the approximately Gaussian zero-mode wavefunctions. Wavefunctions away from 
the Gaussian peak will depend on details of the global structure of X, but the wavefunctions are 
exponentially small in this region. Assuming a Gaussian profile with a fixed width for wavefunc- 
tions is surely a very crude approximation, but it might suffice as a zeroth order approximation. 
A Gaussian landscape on a D = 3-dimensional manifold B assumes that the moduli of stable 
gauge field configurations are so rich that the center coordinates of various fields are scanned 
(almost) randomly and independently. The results of section [S] provide a qualitative picture of 
how the distributions of observables depend on the base manifold B. 

The most important parameter in a Gaussian landscape is the ratio of the Gaussian width 
to the size of the extra dimensions, d/ L. For example, the overall hierarchy of Yukawa couplings 
is proportional to {L/dY on a logarithmic scale. As we have seen in section mi this ratio is 
proportional to Ls/Lg = yo\{B = S'^)/vol(T^). This is generalized to 



As in sections 17.11 and 17.21 the volume of the fiber has to be sufficiently smaller than that of 
the base manifold B in order for hierarchy to be generated. It is an interesting question whether 
this property can be understood within string theory, not as a phenomenological requirement. 

The Heterotic string theory compactified on a T^-fibered geometry corresponds to an eleven- 
dimensional supergravity compactification on a i^3-fibered geometry (this is an intersecting 
D6-D6 system of the Type IIA string theory) in the limit of small vol(T^). Thus, Gaussian 
landscapes on three-folds for various values of (d/L) are intended to simulate a class of vacua 
of string theory that interpolates between the Heterotic theory and M-theory (ignoring stringy 
corrections). The expression in Type IIA language (the last term) in 01831) is known. Numerical 
simulation can tell us the value of d/L that fits to the observed pattern of masses and mixings, 
which suggests the approximate value of (11831) for the vacuum of our universe. 

7.5 T^-Fibered Compactification and F-theory Dual 

The analysis of D = 2 Gaussian landscapes in section [5] may help one understand the landscape 
of Yukawa couplings of Heterotic compactification on a T^-fibered geometry. Let us now consider 
a six- dimensional manifold X that is a T^-fibration over a four-fold B. Let (xg, Xio) parametrize 
a local patch of the T^-fibration and (xs.gj.s) a local patch of B. Fermion zero modes are 
localized on a two-dimensional sub-manifold (called the matter curves) on B so that the gauge 
fields Aq{x5^qj^s) and Aio (2^5,6,7,8) both vanish [31]. Furthermore, their wavefunctions around the 




vol{B) 1 3 1 



(183) 
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matter curves are approximately Gaussian in the two transverse directions determined by VAg 
and VAiQ. The existence of the two transverse directions in which zero- mode wavefunctions are 
Gaussian is quite similar to the D = 2 Gaussian landscapes. One can determine the value of 
d/L that fits the observed hierarchy of masses and mixing angles, which might then be used to 
infer vol(-B)3/vol(T^)5 for the T^-fibered geometry. This class of vacua interpolates Heterotic 
string theory and F(type IIB string)-theory. 

There are also some features that are not captured by the D = 2 Gaussian landscapes. Sup- 
pose that Ag varies along X7, and Aiq along xg] at least it is possible to choose a coordinate 
system on B locally so that this happens. Then a matter curve is along the {x5,Xq) directions. 
However the net chirality of a fermion depends on F^q, the gauge field along the matter curve. 
Therefore zero-mode wavefunctions should have some {x5,Xq) dependence, but no such depen- 
dence is taken into account in the D = 2 Gaussian landscapes. Moreover, the global geometry 
of matter curves may be complicated in B, and they can intersect one another, but the D = 2 
Gaussian landscapes ignore this as well. The behavior of zero-mode wavefunctions on a T^- 
fibered geometry and the geometry of the intersection of matter curves must be studied further 
to determine whether the D = 2 Gaussian landscapes can be useful in studying the landscapes 
of Yukawa couplings with T^-fibered compactification. 

8 Conclusions 

This decade has seen the emergence of a major debate: to what extent is nature fundamentally 
uniquely prescribed, for example by symmetries, vs. to what extent nature results from the 
statistics of a huge landscape of vacua, modified by cosmological and environmental selection. 
For the Standard Model symmetries play a key role, but for physics beyond the Standard Model 
the question remains largely open. Unified gauge symmetries have striking achievements, for 
example a simple interpretation of the quantum numbers of a generation and a precise numer- 
ical prediction for the ratios of the measured gauge couplings. For quark and lepton flavor, 
however, the picture offered by approximate flavor symmetries (AFS) is much less compelling, 
lacking both theoretical simplicity and significant successful predictions. Although the AFS de- 
scription of flavor is apparently well suited to give an understanding of the hierarchical nature of 
charged fermion masses and the CKM mixing matrix, it comes with too much flexibility: with 
an appropriate choice of charges and symmetries any pattern of flavor can be generated. 

In this paper we introduce Gaussian landscapes as some of the simplest landscapes in extra 
dimensions that can account for flavor. Particles are assumed to possess localized (Gaussian) 
zero-mode wavefunctions over some geometry of extra dimensions, and the Gaussian landscape 
consists of independently scanning the peak position of each of these wavefunctions. Small flavor- 
symmetry breaking parameters are replaced by small overlap integrals of these wavefunctions 
on the extra dimensions. Localized zero-mode wavefunctions in extra dimensions are a natural 
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expectation, and may have a more elegant realization than Higgs potentials for flavor symme- 
try breaking. We claim neither precise predictions nor do we present a compelling top-down 
theoretical model, rather we present the patterns of flavor that emerge from these simple toy 
landscapes, the features of these extra-dimensional landscapes that are relevant to the flavor 
problem, and how these landscapes might be realized within the context of string theory. 

In the simplest Gaussian landscape describing the quark sector, where all quarks and the 
Higgs have a universal Gaussian wavefunction (but with independently scanning center coordi- 
nates) over a single extra dimension with geometry 5*^, we flnd the three major characteristics 
of quark flavor: a hierarchical distribution of quark masses, pairing structure (the VT-current 
approximately connects distinct pairs of quarks), and generation structure (the electroweak 
pairing connects the heaviest up-type quark to the heaviest down-type quark, and similarly for 
the middle and lightest quarks). The relevant probability distributions are shown in Figures H] 
and [5] and result from inputting just two parameters: a universal constant g^s of order unity 
setting the scale for overlap integrals and a universal constant d/L setting the width of Gaussian 
wavefunctions relative the size of the extra dimension. 

The flavor structure of the quark sector can also be obtained from AFS using two free 
parameters. For example, by hypothesizing an approximate U{1) symmetry with a leading 
Yukawa coupling of order unity and others suppressed by various powers of a small symmetry 
breaking parameter e. However there is a crucial difference. In the AFS case one must carefully 
choose the U{1) charges of each of the flfteen fermions of the Standard Model. A huge variety of 
mass patterns could be accommodated by suitable charge choices. In the Gaussian landscape no 
such choices are made. Each of the flfteen fermions is treated symmetrically. They differ only 
by the location of their Gaussian wavefunction, and these are scanned randomly over S^. Thus 
the hierarchies arise purely from statistics; they cannot be changed as they do not involve any 
free parameters beyond Qcs and d/L. 

While the above accomplishments of the Gaussian landscape on are striking, there are 
certain features that are less than ideal. Although they are peaked, the probability distributions 
for flavor parameters are quite broad, as can be seen in Figure [5l At half maximum, the 
deviation from the peak value is typically an order of magnitude. Thus the statistical nature of 
the landscape prevents us from making precise predictions. We flnd that this order of magnitude 
width is also typical of Gaussian landscapes in more than one dimension. The landscape makes 
no distinction between up and down sectors. This is a problem for explaining the observed large 
t/b mass ratio, in particular because on the distribution for the top and bottom masses is 
narrower than for the other generations and also because it is peaked near the maximum value. 
Thus even if there were a selection effect favoring a heavy top quark, the observed bottom mass 
would still be somewhat unlikely. 

One possibility to resolve this shortcoming is to replace the effective coupling g^s with two 
parameters: one for the up sector that is about an order of magnitude larger than the one for 
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the down sector. It would be interesting to find the landscape origin for such a "tan factor. 
Another possibility is that the shape of the t/h mass distribution is a special feature of the 
landscape. We find that this is the case. In the landscape, where Gaussian wavefunctions are 
distributed at random over the surface of a (square) torus, the distribution is wider, as shown in 
Figure [TTl The same result holds for Gaussian landscapes on S"^ and . A selection for a heavy 
quark, for example for electroweak symmetry breaking, could then more easily account for the 
t/h mass ratio. For any number of extra dimensions, we can analytically compute the Yukawa 
probability distribution near maximal values, and we find that it is suppressed as the dimension 
increases, strengthening this interpretation of the t/h ratio. 

Although flavor symmetries play no fundamental role in Gaussian landscapes, small overlap 
integrals may allow AFS to emerge in the low-energy theory. For example, if gi, q2 and h have 
narrow Gaussian profiles with centers that are well-separated from each other, then AFS that 
act on qi and q2 will emerge at low energy. We have found that Yukawa matrices have the form 
expected from Abelian AFS, namely \ij ~ e^e^, if the Higgs width is not too wide and if the 
Yukawa coupling is not too small. As the separation between the Higgs and quarks approaches 
the scale of the extra dimension, effects associated with the periodicity of the wavefunctions 
destroys the AFS form of the Yukawa coupling. Nevertheless, the AFS approximation aids 
in understanding the numerical results, and we have used it to compute approximate analytic 
distributions for the AFS factors e, the quark mass eigenvalues, and the CKM mixing angles. 

This strategy is particularly useful for understanding Gaussian landscapes on higher dimen- 
sional geometries where numerical integration requires greater computational resources. Analytic 
results for the geometries S'^, S'^, = S*^ x 5*^, and have been obtained and the distributions 
are compared in Figure [T5l We find that the qualitative predictions of the Gaussian landscape on 

remain intact. One important feature is that, independent of the dimension D, the AFS dis- 
tribution functions and flavor observables are found to be polynomial or logarithmic functions 
of InAj or In(sin^jj). By comparing numerical distributions of flavor parameters in Gaussian 
landscapes on and 5^, we have investigated the importance of the shape and curvature of 
the extra dimensional manifold. We find some differences in the details, but the more striking 
are the similarities, as seen for example in Figure [TH 

Large leptonic mixing angles suggest that in the low-energy theory no AFS emerges for the 
lepton doublets Zj, implying that the Gaussian widths for U are comparable to the size of the 
extra dimension(s). Since the fermion mass hierarchies are smaller in the down and charged 
lepton sectors than in the up sector, we economize on parameters and assume that the widths 
preserve an SU{5) symmetry, with ^ dio, du. The smallness of the neutrino masses suggests 
that there is no light right-handed neutrinos, such that the light neutrinos are Majorana and the 
origin of neutrino masses involves the breaking oi B — L symmetry. In the seesaw mechanism 
this implies Majorana masses for heavy right-handed neutrinos. In the Gaussian landscape we 
expect Gaussian profiles for both Vji and the B — L breaking fields. The center of the B — L 
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breaking is not expected to be correlated with the center of SU{2) x U{1) breaking, leading to 
a very different statistical character for the Majorana and Dirac mass matrices of neutrinos. 

We find that the physics of neutrino masses and mixings is quite different depending on 
whether the profile for 17 fi is narrow or wide. In the case that it is narrow, localization leads to 
significant hierarchies in the eigenvalues of both the Majorana and Dirac matrices. In theories 
based on AFS these hierarchies cancel in the light neutrino mass matrix, but in the Gaussian 
landscape these hierarchies add. The result is that 7712 /m^ is typically too small to agree with 
observation, unless perhaps if strong selection effects are important. For wide V^i profiles, the 
neutrino Yukawa and Majorana matrices take on a democratic form, and again have hierarchical 
eigenvalues. If the profiles are sufficiently wide, 7712/ is always less than the observed value, 
as shown in Figure [T^ Similar difficulties arise with the lepton mixing angle 623- Although a 
large width for the lepton doublets leads to large 23 rotations to diagonalize the charged and 
neutral lepton mass matrices, these large angles cancel so that the physical mixing angle 6*23 is 
typically small, as shown in the numerical simulation of Figure [21] 

These difficulties, that 777,2/^3 and 6'23 are typically much smaller than the values observed, 
are both solved by a remarkably simple observation. In order to account for CP violation in 
the quark sector a complex phase must be introduced to the Gaussian profile. Using the simple 
profile of flSSl) . we find the following results for the Gaussian landscape on S^: 

• Previous results for distributions of charged fermion masses are preserved, for example 
Figure [16] is essentially unchanged for CP- violating Gaussian landscapes. 

• The generation structure of the quark sector is preserved (Figure [2^ . 

• The CKM phase is of order unity (Figure [25]) . 

• Large leptonic mixing angles 6*12 and 6'23 are obtained (Figure [23]) . 

• The large observed value for the neutrino mass ratio 7712/ is not atypical (Figure [^UT h)). 

The distributions of flavor observables on Gaussian landscapes are typically broad. How- 
ever, given the experimental measurement of some subset of flavor observables, the conditional 
distributions for the remaining observables in our universe changes significantly. We find that 
even quite loose cuts on the measured neutrino parameters are sufficient to considerably sharpen 
the predictions for the leptonic mixing angle, ^13, and the neutrinoless double beta parameter, 
m^^, within a particular Gaussian landscape; see Figure [2BI and fllU2p . respectively. Furthermore 
these cuts lead to a prediction of large CP violation in neutrino oscillation as shown in the last 
row of Figure [25] and in the scatter plots of Figure [26] 

When a supersymmetric gauge theory is compactified on an internal manifold, quarks, leptons 
and the Higgs boson may originate from gauginos and the components of gauge fields with 
polarization along the internal manifold. Their Yukawa interactions originate from the gauge 
interactions of the gauginos, while the Yukawa coupling constants are calculated by overlap 
integration of zero-mode wavefunctions over the extra-dimensional space. Wavefunctions become 
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Gaussian when the extra-dimensional manifold is a torus fibration, and the fiber is small relative 
to the base. Thus, Gaussian landscapes fit very well into this framework — compactification of a 
supersymmetric gauge theory on a torus-fibered manifold. 

The ratio d/L, which controls the hierarchy, roughly reflects the ratio of the torus fiber size to 
the base size in Heterotic theory language. We now know the value oi d/ L that is good for quark 
and lepton phenomenology, but it remains an open question why such a value has been selected 
from all the possibilities in the landscape (or if it represents a typical subset of the landscape). 
Indeed, there is no theoretical motivation to consider torus-fibered compactification over the 
myriad of other geometries for compactification. We also assumed in the Gaussian landscape 
that the center coordinates of all Gaussian wavefunctions scan randomly and independently over 
the base space. This assumption is not totally without motivation, but it would be nice if the 
validity of this assumption were studied in string theory. 

Simple Gaussian landscapes, motivated by supersymmetric gauge theories compactified on 
a torus-fibered manifold, can account for the broad pattern of quark and lepton masses and 
mixings in terms of just five parameters. While the probability distributions for the measured 
parameters are broad, predictions from a given landscape for future measurements in neutrino 
physics are more precise. 
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A Approximate Probability Distribution Functions 

We here collect the probability distribution functions calculated using the AFS approximation, 
i.e. using and (13BHIU]) . for the Gaussian landscape on S^. First we list the distribution 
functions of the Yukawa eigenvalues that follow from fl33p and the approximation fl36p : 




for < Z3 < 1 



dPiz^ 

dZ3 



-— 32 - 76v/i^^ + 48^3-26z3v/i^^ 



= < 



(184) 



- 3 (8 + 8^3 + ^1) arccot {^/zi^) 



- 3 {8 + 8z3 + z^) arctan [y/z^ - l) 



for 1 < Z3 < 2 . 
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r 3 



dP{z2) 

dZ2 



Z2 (24 - + Sttz: 



for < Z2 < 1 



- (4 - 2y^I^ - Qz2 + 5^2^^^^) 



;i85) 



dPjzi] 
dzi 



9 2 

4 ^ 



arctan (V^2 — l) — arccot (V 22 — l) 



^|2(2-^i)v/ir^ 

arctan (a/^i — l) — arccot — l) 



for 1 < ^2 < 2 . 

for < 21 < 1 , 

for 1 < zi < 2 . 



(186) 



Here zt = ln(A,/Amax)/Alne where A^ax = (4LV97rrf2)i/4^^g and Alne = -{L/d^/U. We also 
note the mean values of these distributions: 



(Z3) = 0.2, (^2) = 0.6, (^i) = 1.2. 



(187) 



Meanwhile, the distribution functions of the CKM mixing angles follow from (133|1 and the 
set of approximations (ISBHIU]) . They are given by: 



dP{t 



12 



dti2 

dP{t23) 



6(l-v^) (l + 2v^)c?ti2 



dt 



23 



-2^1 — ^23 + arccosh ( 1 

^23 



X 



2VI - ^23(1 + 2t23) - 3t23arccosh 



^23 



dP{t 



13) 



dt 



13 



6(v^r^tl^+v^-l) (1- v^) 



X 



2ti3 + 2v/ti3(l - t 



13 j 



(189) 



(190) 



where each applies over the full range < tij < 1, with tij = In | sin 6'jj| /Alne. The mean values 
of these distributions are: 



(ti2) = 0.17, (t23) = 0.10, (ti3)=0.35. 
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